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ABSTRACT 



Background: Transcranial direct current stimulation (tDCS) is a form of neuromodulation that is in- 
creasingly being utilized to examine and modify a number of cognitive and behavioral measures. The 
theoretical mechanisms by which tDCS generates these changes are predicated upon a rather large 
neurophysiological literature. However, a robust systematic review of this neurophysiological data has 
not yet been undertaken. 

Methods: tDCS data in healthy adults (18-50) from every neurophysiological outcome measure reported 
by at least two different research groups in the literature was collected. When possible, data was pooled 
and quantitatively analyzed to assess significance. When pooling was not possible, data was qualitatively 
compared to assess reliability. 

Results: Of the 30 neurophysiological outcome measures reported by at least two different research 
groups, tDCS was found to have a reliable effect on only one: MEP amplitude. Interestingly, the mag- 
nitude of this effect has been significantly decreasing over the last 14 years. 

Conclusion: Our systematic review does not support the idea that tDCS has a reliable neurophysiological 
effect beyond MEP amplitude modulation - though important limitations of this review (and conclusion) 
are discussed. This work raises questions concerning the mechanistic foundations and general efficacy of 
this device - the implications of which extend to the steadily increasing tDCS psychological literature. 

© 2014 Published by Elsevier Ltd. 
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1. Introduction 

Nearly fifteen years ago, researchers revived a field of research 
that explores the effects of running a weak electric current be- 
tween two electrodes placed on the scalp of healthy individuals 
(Konig and AnkenmuUer, 1960; Hamer, 1968; Lolas 1977; Elbert 
et al., 1981; Nitsche and Paulus, 2000). Dubbing this technique 
transcranial direct current stimulation (tDCS), experiments 
showed a marked, time dependent, and polarity specific modula- 
tion of neuronal populations underlying the electrodes. Since this 
initial neurophysiologic finding, a growing number of researchers 
and clinicians have been exploring the effects of tDCS across a 
number of cognitive/behavioral domains. Today, the literature 
suggest tDCS can enhance a number of higher-order cognitions 
and behaviors ranging from working memory and motor learning 
to emotional regulation and focused attention (for review: Nitsche 
and Paulus, 2011). The claims made by cognitive and behavioral 
tDCS researchers largely depend upon the mechanistic framework 
suggested by the neurophysiologic data. However, a robust sys- 
tematic review of the neurophysiologic impact of tDCS has not yet 
been undertaken. This is something we hope to remedy in this 
paper. 

1.1. tDCS: A brief overview and proposed mechanisms of action 

Modern tDCS devices typically consist of an adjustable direct 
current stimulator and two stimulating electrodes (an anode and a 
cathode). These electrodes are typically attached to two separate 
locations on the scalp (either directly or via larger sponge elec- 
trodes) and a weak current (0.5-2.0 mA) is run between the 
electrodes. As this current passes between the electrodes, it is 
believed a small amount of the current passes through the brain. 
This current flow is purported to modulate neural activity 



underneath the electrode and, to a lesser extent, diffuse locations 
in the brain (Nitsche et al., 2008; for debate: Bikson. 2013). 

There are two mechanisms by which tDCS modulates brain 
activity that are widely accepted in the field. The first proposes 
tDCS modulates the resting membrane potential of neuronal po- 
pulations via ionic adjustment of extracellular space. More speci- 
fically, neurons proximal to the anode are thought to become 
hypo-polarized whist neurons near the cathode are thought to 
become hyper-polarized (Stagg and Nitsche, 2011). This shift in 
resting membrane potential is believed to occur both during sti- 
mulation and for a short period of time ( < 5 min) following sti- 
mulation. The second proposes tDCS modulates synaptic activity in 
a manner akin to long term potentiation (under the anode) and 
long term depression (under the cathode: Stagg and Nitsche, 
2011). This mechanism is believed to be active for an extended 
period of time (up to 120 min) following the cessation of long- 
duration ( > 7 min) stimulation. In this systematic review, we 
group studies into short- and long-duration stimulation to account 
for membrane and synaptic effects believed to be triggered by 
tDCS. 

1.2. Systematic review structure 

The neurophysiological effects of tDCS have largely been mea- 
sured utilizing four approaches: transcranial magnetic stimulation 
(TMS), event related potentials (ERPs), electroencephalographic 
spectral analyses (EEC), and functional magnetic resonance ima- 
ging (fMRl). Accordingly, the methods and results sections (below) 
will be structured around these modalities. Additionally, the ma- 
jority of studies exploring the neurophysiological effects of tDCS 
have utilized a single measure: TMS motor evoked potential (MEP) 
amplitude. Due to the disproportionately large number of studies 
exploring this measure, we have decided to dedicate an analysis 
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1 section solely to this measure (rather than combine it with other refractoriness and GABAergic mediated synaptic inhibition (Na- 67 

2 TMS measures). l<amura et al., 1997). When the two pulses are applied with an ISI 68 

3 of 7-30 ms, the second MEP will typically be larger than the first. 69 

4 1.2.1. Transcranial magnetic stimulation (TMS) This elevation, termed ICF, is thought to reflect glutamatergic 70 

5 Transcranial magnetic stimulation (TMS) is a device that allows mediated synaptic facilitation. During SICl and ICF paradigms, test 71 

6 for non-invasive electrical stimulation of localized regions of the pulse MEP amplitudes are typically reported as a percentage of the 72 

7 cortical surface. Via the rapid passing of a strong electrical current normalized conditioning pulse MEP amplitudes (Nakamura et al., 73 

8 through a coil held against an individual's head, TMS generates an 1997). 74 

9 equally rapid fluctuating magnetic field which (via electro- As each of these TMS outcome measures have been utilized to 75 

10 magnetic induction) generates a secondary current within the explore the modulatoiy effects of tDCS, the results of each will be 76 

11 cortex effectively force-firing underlying neuronal populations. presented in turn. 77 

12 Single TMS pulses of sufficient strength to the motor cortical re- 78 

13 presentation of a particular muscle can cause that muscle to ac- 1.2.2. £vent related potential (ERP) 79 

14 tivate. Various aspects of muscular activation can be measured ERPs are specific and reliable neural responses to specific sen- 80 

15 (using electromyography: EMG) and analyzed to infer cortical and sory stimuli. These responses can be measured using electro- 81 

16 spinal excitability patterns (Hallett, 2007). As such, TMS evoked encephalography (EEC). A typical ERP will contain a number of 82 

17 muscular responses represent an excellent tool to explore mod- negative and positive peaks occurring within predictable time 83 

18 ulation engendered by tDCS. There are several TMS evoked motor windows following stimulus onset. The constituent parts of ERPs 84 

19 outcomes which have been utilized to explore the modulatoiy are thought to reflect brain function. More specifically, peak am- 85 

20 effects of tDCS. plitude is thought to reflect underlying cortical excitability (Handy, 86 

21 2005). As such, the comparison of ERP amplitudes prior to and 87 

22 1.2.1.1. Motor evoked potential (MEP) amplitude. When a TMS pulse following tDCS represents a unique way to measure the mod- 88 

23 activates a target muscle, the resultant muscular 'twitch' is com- ulatoiy effects of tDCS. Although a number of ERP measures have 89 

24 monly referred to as a motor evoked potential (MEP). The ampli- been explored in the literature, only four have satisfied our 90 

25 tude of each MEP can be thought of as a measure of the strength of inclusion criteria of being replicated by at least two different 91 

26 the muscular activation. As it is generally accepted that TMS ac- research groups (for a full list of ERP measures that did not meet 92 

27 tivates pyramidal neurons transynaptically (Rotenberg et al., 2014), this criteria, see Table SI ). 93 

28 TMS generated MEP amplitudes are believed to reflect the excit- 94 

29 ability of corticocortical circuits. Additionally, MEP amplitudes are 1.2.2.1. PlOO visual evoked potential (VEP). When an ERP is evoked 95 

30 reflective of and can be modulated according to, spinal mo- using visual stimuli (either low or high contrast), the resultant 96 

31 torneuron excitability. Accordingly, when utilizing a single TMS response is termed a VEP. VEPs often include a positive peak 97 

32 pulse intensity, any increase in MEP amplitude is thought to reflect ~ 100 ms post stimulus presentation at the occipital pole. Termed 98 

33 either an increase in synaptic excitability at the level of the motor the PlOO, the amplitudes of this peak has been shown to be highly 99 

34 cortex or the spinal cord- and vice versa (Hallett, 2000). reliable and sensitive to excitability shifts within individuals at the 100 

35 level of the visual cortex (Di Russo et al., 2002). QSOl 

36 1.2.1.2. Resting motor threshold (rMT) and active motor threshold 102 

37 (aMT). rMT commonly refers to the minimum intensity of TMS 1.2.2.2. N20 somatosensory evoked potential (SEP). When an ERP is 103 

38 pulse required to elicit a specified number of MEPs (typically 50% evoked using a brief sensory stimulus (typically generated on the 104 

39 response rate) in a relaxed muscle. Although synaptic connections peripheiy of the body), the resultant response is termed an SEP. An 105 

40 at the cortical level doubtless influence this measure, rMT is often utilized method of SEP evocation is the use of electrical 106 

41 thought to depend more on axonal membrane excitability (Hallett, stimulation at the median nerve of the wrist. This particular SEP 107 

42 2007). aMT commonly refers to the minimum intensity of TMS often includes a negative peak ~20 ms post-stimulation near the 108 

43 pulse required to elicit a specified number of MEPs (typically 50% somatosensoiy cortex. Termed the N20, the amplitude of this peak 109 

44 response rate) in an activated muscle (typically voluntary con- has been shown to be highly reliable and sensitive to excitability 110 

45 traction at 15-30%). As muscle contraction elevates the excitability shifts within individuals at the level of the somatosensory cortex 111 

46 of large groups of cortical neurons and spinal motoneurons, aMT is (Emerson et al., 1988). 112 

47 thought to depend more directly on axonal threshold than rMT 113 

48 (Hallett, 2007). 1.2.2.3. N2 and P2 laser evoked potential (LEP). When an ERP is 114 

49 evoked using a heat generating laser (typically focused on the back 115 

50 1.2.1.3. Cortical silent period (cSP). ATMS pulse to the motor cortex of the hand), the resultant response is termed an LEP. Thought to 116 

51 typically engenders both an excitatory and inhibitoiy effect. represent neural pain processing, LEPs often include a negative/ 117 

52 Whereas the excitatory effect is reflected in the MEP, the inhibitory positive wave ~200 ms post-stimulation near the somatosensory 118 

53 effect is reflected in a suppression of EMG activity ( ~ 100-300 ms) cortex. Termed the N2 and P2 waves (respectively), the amplitudes 119 

54 immediately following the MEP; the duration of which constitutes of these peaks have been shown to be highly reliable and sensitive 120 

55 the cSP. The first 50 ms of the cSP are thought to reflect spinal to excitability shifts in individuals at the level of the somatosen- 121 

56 mechanisms, whereas the remaining portion is mediated by cor- sory cortex (Ohara et al., 2004). 122 

57 tical inhibition (Cantello et al., 1992). 123 

58 1.2.2.4. - Mu/Alpha event related desynchronization (ERD). 124 

59 1.2.1.4. Short interval cortical inhibition (SICI) and intracortical fa- Although not a singular potential, ERD is a measure of amplitude 125 

60 cilitation (ICF). The application of two TMS pulses to the same response across an oscillatory range (see below for additional EEG 126 

61 motor cortical location separated by variable interstimulus inter- power spectrum information). Amplitude modulation across un- 127 

62 vals (ISIs) is another method for exploring excitatory and in- ique frequencies has been shown to correlate with unique cogni- 128 

63 hibitory corticocortical circuits. When the first (conditioning sti- tive and/or behavioral activity. More specifically, the Mu/Alpha 129 

64 mulus) and second (test stimulus) TMS pulses are applied with an frequency (typically 8-13 Hz) over the sensoiy-motor cortex ty- 130 

65 ISI of 2-6 ms, the second MEP will typically be smaller than the pically displays large amplitude peaks when a person is at rest. 131 

66 first. This diminution, termed SICI, is thought to reflect axonal However, these waves attenuate (de-synchronize) during motor 132 



Please cite this article as: Horvath, J.C, et al.. Evidence that transcranial direct current stimulation (tDCS) generates little-to-no reliable 
neurophysiologic effect beyond MEP amplitude.... Neuropsychologia (2014), htlp://dx. doi. org/10. 1016/j. neuropsychologia. 2014.1 1.021 



ARTICLE IN PRESS 



4 J.C. Horvath et al. / Neuropsychologia I fllllj ■■■-■11 

1 performance, obsei"vation, and imageiy - theoretically due to overt significant effect of tDCS on this measure whilst two reported no 67 

2 or covert activation of the motor system (Guieu et al., 1999). Al- effect. 68 

3 though this outcome measure may be more comparable to those Two papers have explored the impact of tDCS on bioenergetics 69 

4 included in the EEC power spectrum section (below), we have using Phosphorus MRS (Binkofski et al., 2011; Rae et al., 2013). 70 

5 included it here as we were able to analyze it quantitatively. Unfortunately, the stimulation protocol and outcome measures 71 

6 reported for each were different negating the possibility of a 72 

7 3.2.3. Electroencephalography (EEC) power spectmm meaningful comparison. 73 

8 EEC is an imaging modality used to record the electrical activity Two papers have explored the effects of tDCS on regional cer- 74 

9 of neural groups. Although EEC recordings can be divided into and ebral blood flow using positron emission tomography (Lang et al., 75 

10 analyzed according to brief temporal windows (see ERPs, above), 2005; Paquette et al., 2011). Again, as the tDCS protocols utilized 76 

11 long duration EEC analysis allows for the determination of oscil- by each different significantly, it was not possible to generate a 77 

12 lation patterns over-time across the entire brain or within sped- meaningful comparison. 78 

13 fied neural regions. Oscillatory patterns are typically broken into Finally, several recent studies have utilized HD-tDCS and in- 79 

14 five different frequency bands; Delta (0-4 Hz), Theta (4-7 Hz), dividualized electrode placement based on personal MRI electric 80 

15 Alpha (7-14 Hz), Beta (15-30 Hz), and Gamma ( > 30 Hz). These current flow models (e.g. - Kuo et al., 2013; Villamar et al., 2013; 81 

16 oscillations represent synchronous activity over diffuse neural Borckardt et al., 2012; Clark et al., 2012). Unfortunately, although 82 

17 networks and the power of each (relative to the others) is thought these new techniques are promising, most studies have utilized 83 

18 to correspond to different states of brain function and/or behavior clinical population and there is not yet enough data in healthy 84 

19 (e.g.: high power in the theta band has been associated with controls to derive meaningful conclusions. Accordingly, we have 85 

20 drowsiness and response inhibition: Rowan and Tolunsl<y, 2003). omitted these studies from our analysis. 86 

21 As such, the unique ratio of oscillatoiy power prior to and fol- 87 

22 lowing stimulation whilst individuals undertake varied tasks re- 88 

23 presents an important method for exploring the modulatory effect 2. Methods 89 

24 of toes on diffuse cortical networks. 90 

25 2.1. General 91 

26 1.2.4. Functional magnetic resonance imaging (fMRI) 92 

27 fMRI is used to measure blood oxygen level dependent (BOLD) It is common in quantitative analyses of multiple studies to 93 

28 signals within the brain across time. The hemodynamics measured translate the outcome values of each to a common measure before 94 

29 by fMRI are thought to correspond with neural activity such that undertaking any comparison (e.g. - effect size or standardized 95 

30 an increase of blood flow to a cortical region is thought to be as- mean difference). In order to obtain a common measure, however, 96 

31 sociated with an increase in neural activity within that region reported outcome values must be compared to a control or sham 97 

32 (Logothetis et al., 2011). The amount or intensity of BOLD activa- condition. Unfortunately, the vast majority of numerically analyz- 98 

33 tion prior to and following stimulation can be used to explore the able studies included in this systematic review (92 out of 117) did 99 

34 resting firing rate modulatory effects of tDCS both underneath the not include a control or sham condition. For this reason, no effect 100 

35 electrode and across diffuse cortical networks. . size and/or standard mean difference could be determined. Rather, 101 

36 we maintained the central measure as reported in each paper and 102 

37 1.3. Omitted measures _ ■■^0'"^ obtained weighted means and pooled standard deviations for 103 

38 normalized data at each appropriate time-point post-stimulation. 104 

39 For a list of TMS, ERP, EEG power spectrum, and fMRI related If numerical data was not supplied in the text, averages and var- 105 

40 outcome measures not replicated between at least two research lance were extracted from included images (achieved by exporting 106 

41 groups (therefore, excluded from this analysis), please see Table SI images to an image editing program, overlaying a standardized 107 

42 in the Supplementary material. Several additional physiological grid, and counting relevant values by hand). This visual data ex- 108 

43 outcome measures (not easily grouped into one of these 4 cate- traction process was undertaken by the lead author. When possi- 109 

44 gories) were also not included in this analysis due to using non- ble, obtained values were run through the statistical tests reported 110 

45 comparable protocols and/or data reporting. Additional protocol in the representative paper to assess accuracy of the extracted 111 

46 details from each of these studies (briefly outlined below) are data. 112 

47 presented in Table S2 in the Supplementary material. Final values were compared using unpaired Student's f-tests 113 

48 Although a number of papers have used MRI to explore func- (when no SHAM condition was available), and one-way ANOVAs 114 

49 tional connectivity measures following tDCS to Ml (Polania et al., (when a SHAM condition was available) to determine any differ- 115 

50 2012a, 2012b; Sehm et al., 2013; Stagg et al., 2009b) and the DLPFC ence between obtained values. If significance using an ANOVA was 116 

51 (Keeser et al., 2011; Park et al., 2013; Pena-Gomez et al., 2012), the established, post-hoc analysis consisted of between condition 117 

52 modes of data acquisition, analyses, and reporting conducted by bonferonni corrected f-tests. As each experiment reported varied 118 

53 each are dissimilar and difficult to meaningfully compare. For temporal testing protocols following stimulation (see tables for 119 

54 these reasons, we have opted not to include these measures in this more specific information), we were unable to run any two-way 120 

55 systematic review. It is worth noting, however, that Paquette et al. ANOVAs (using 'time' as a variable). 121 

56 (2011) recently noted that these papers have produced "...to some As noted below, one of the only absolute inclusion criteria for 122 

57 extent, conflicting results" (p. 2086). this analysis was that an outcome measure be replicated by at 123 

58 Additionally, four papers have used magnetic resonance spec- least two different research groups. Examples of non-experimen- 124 

59 troscopy (MRS) to explore the neurochemical impact of tDCS tally based systematic errors reliably influencing the results gen- 125 

60 (Clark et al., 2011; Rango et al., 2008; Stagg et al, 2011, 2009a). erated by a single research group can be found in all the sciences, 126 

61 Unfortunately, as each used different and incompatible imaging from physics (see Gieryn, 1992), to biology (see Mobley et al., 127 

62 methodologies and stimulation parameters, we have opted not to 2013), to medicine (see Osherovich, 2013). We, therefore, chose to 128 

63 include these measures in our systematic. It is worth noting, exclude measures that have only been replicated by a single re- 129 

64 however, that the only common neurometabolite measure re- search group to ensure all data included in and conclusions gen- 130 

65 ported between the studies was Glx concentration (a composite erated by this review accurately reflect the effects of tDCS itself 131 

66 measure of glutamate and glutamine). Of the four, two reported a rather than any unique device, protocol, or condition utilized in a 132 
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Table 1 

Studies exploring tiie effects of low and higli density tDCS on TMS elicited MEP amplitude of varied muscles at rest. 



Paper/* of studies 


N 


tDCS duration 


tDCS polarities 


DB/NN 


Tar. 


TMS power/frequency 


# of MEPs averaged 


Tlx) 


Low density (0.02857 mA/cm^)/very-short duration ( < 1 min) — tested muscle at rest 












Nitsche and Paulus (2000) 


10 


4 s 


A, C 


-/- 


IHM 


2 mV MEP/0.25 Hz 


12 


- 


Nitsche et al. (2003a, 2003b) 


12 


4s 


A, C 


-/- 


IHM 


1 mV MEP/0.1 Hz 


15 


- 


Nitsche et al. (2004) 


6 


4s 


A, C 


-/- 


IHM 


1 mV MEP/0.1 Hz 


15 


- 


Nitsche et al. (2004) 


12 


4 s 


A. C 


-/- 


IHM 


1 mV MEP/0.1 Hz 


15 


- 


Nitsche et al. (2007/x3) 


12 (x3) 


4 s (x3) 


A, C (x3) 




IHM 


1 mV MEP/0.1 Hz 


15 




Low density (0.02857 mA/cm^)/short duration (1- 


7 min) — tested muscle at rest 














Nitsche and Paulus (2000) 


19 


7 min 


A, C 


-/- 


IHM 


2 mV MEP/0.25 Hz 


Continuous 


TO-TIO 


Nitsche and Paulus (2001 ) 


12 




/\ 


-/- 


IHM 


1 mV MEP/0.25 Hz 


20 bl/15 @ 5 min 


T0-T15 


Liebetanz et al. (2002) 


11 


5 min 


A. C 


Y/- 


IHM 


1 mV MEP/0.25 Hz 


20 bl/continuous 


TO-TIO 


Nitsche et al. (2003/x2) 


12 (x2) 


5 min (xl)7 min (xl) 


C (x2) 




IHM 


1 mV MEP/0.25 Hz 


20 bl/15 @ 5 min 


TO (xl)TO-TlO (xl) 


Nitsche et aL (2004) 


6 


7 min 


A. C 


-/- 


IHM 


1 mV MEP/0.25 Hz 


20 bl/15 per 5 min 


T0-T15 


Nitsche et aL (2004) 


9 


5 min 


A, C 


-/- 


IHM 


1 mV MEP/0.25 Hz 


20 bl/continuous 


TO-TIO 


Quartarone et al. (2004/x2) 


7(xl) 
21 (xl) 


5 min (x2) 


A (xl) C (xl) 


-/- 


IHM 


0.5-1 mV MEP/0.25 Hz 


15 @ 10 min 


TO, TIO (x2) 


Nitsche et aL (2007/x3) 


12 (x3) 


7 min (x3) 


A. C (x3) 


-/- 


IHM 


1 mV MEP/0.25 Hz 


15 @ 5 min 


T0-T15 


Nitsche et aL (2007) 


12 


7 min 


A, C 


-/- 


IHM 


1 mV MEP/0.25 Hz 


20 bl/15 @ 5 min 


T0-T15 


Fricke et al. (2011/x2) 


8(xl) 
12 (xl) 


5 min (xl)7 min (xl) 


A (x2) 


-/- 


IHM 


1 mV MEP/0.25 Hz 


30 @ 5 min 


T0-T5 (xl)T0-T15 (xl) 


Munneke et aL (2011) 


10 


7 min 


C 




IHM 


0.5 mV MEP/0.25 Hz 


30 


T5 


Schade et al. (2012/x2) 


8 (x2) 


5 min (x2) 


A. C (x2) 


Y/- 


IHM 


1 mV MEP/0.25 Hz 
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Rizzo et al. (2014) 
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Table 1 (continued ) 



Paper/* of studies 



tDCS duration 



tDCS polarities 



DB/NN Tar. TMS power/frequency # of MEPs averaged 



Tlx) 



O 3 

3 

3 



Hasan et al. (2012) 
Suzulti et al. (2012) 
Monte-Silva et al. (2013/x2) 
Hasan et al. (2013) 
Batsikadze et al. (2013/x2) 

Schabrun et al. (2013/x3) 



Simis et al. (2013) 
Teo et al. (2014) 

High density ( > 0.02857 mA/cm^)/long 
Ardolino, Bossi et al. (2005) 

Jeffery et al. (2007/x2) 
Furubayashi et al. (2008) 
Jefferson et al. (2009/x3) 

Groppa et al. (2010) 
Roche et al. (2011) 
Kuo et al. (2013) 
Batsikadze et al. (2013) 
Galvez et al. (2013/x2) 
Tremblay et al. (2013) 
Schmidt et al. (2013) 
Pellicciari et al. (2013) 
Miyaguchi et al. (2013) 
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General notes: DB= double-blind experimental design with regards to tDCS condition (not any included pharmacological arm); NN=neuronavigation system used to guide TMS; Tar.=target muscle tested; T(x)=time (in minutes) 
post-stimulation that measurements were reported and pooled for this analysis. A=anode; C=cathode; S = sham; IHM^intrinsic hand muscle (first dorsal interosseous; abductor digiti minimi; abductor pollicis brevis) LM=leg 
muscle (Tibialis Anterior; Vastus Lateralis); PM=pharyngeal muscle; 1NF= infraspinatus; BB=bicep brachii; ECR=extensor carpi radialis; SO=stimulator output; rMT=resting motor threshold; aMT=active motor threshold. 
TMS power refers to the TMS output utilized in each experiment. For example, 1 mV means the stimulator was set to a level which consistently evoked peak-to-peak MEP amplitudes of 1 mV. 

Low intensity short-duration notes: two studies from Nitsche and Paulus (2000) and one study from Nitsche and Paulus (2001 ) were excluded due to no SD or SEM being discernable from reported data. One study from Fricke et al. 
(2011) was excluded due to SD over 100%. A sham condition from Quartarone et al. (2004) was excluded due to reporting MEP modulation during motor imagery only. 

Low intensify long-duration notes: One study from Nitsche and Paulus (2001) was excluded due to no SD or SEM being discernable from reported data. Boros et al. (2008) and Cambieri et al. (2012) were excluded due to SD over 
100%. McCambridge et al. (2011) was excluded due to no timing of post-stimulation MEP measurements being discernable from reported data. One study from Monte-Silva et al. (2009) was excluded due to it including data 
reported in a prior report. Two studies from Antal et al. (2007) were omitted due to including a cognitive and/or motor task during tDCS stimulation. Four studies from Miyaguci et al. (2013) were omitted due to including a passive 
or active motor activity during stimulation. One study from Bradnam et al. (2010) was omitted due to targeting a muscle under voluntary isometric contraction. 

High intensity notes: one study from Jefferey et al. (2007) was excluded due to targeting a muscle under voluntai'y isometric contraction. Paquette et al. (2011) and Kim et al. (2013) were excluded due to SD over 100%. 
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single lab. Due to this inclusion criterion, a number of outcome 
measures extant in the literature were omitted due to only being 
reported by a single research group (see Table SI). We want to 
emphasize that this does not suggest said research is in any way 
faulty or incorrect; rather, that inter-group replication as necessary 
in order to eliminate any potential non-tDCS influential outcome 
factors. 

Additional information for each measure is outlined below. 
2.2. TMS: MEP amplitude 

Papers included in this analysis were obtained from a PubMed 
database search (January 20, 2014). The search term "transcranial 
direct current stimulation" generated 1040 papers. Initially, the 
abstract of each paper was read to determine outcome measures. If 
a TMS measure was included, the methods section was read to 
determine what protocol/s were utilized. 

The only absolute inclusion criteria for this measure were MEPS 
must have been obtained in healthy adults (18-50), the muscular 
representation of the tested muscle must have been covered with 
the active electrode (this excludes measure of ipsilateral muscular 
effect), the tested muscle must have been at rest (not voluntarily 
contracted) and the reported SD must have been below 100% (this 
final quality-assessment criterion was included to ensure ex- 
cessively variable results did not inequitably influence pooling). 
Any study which fulfilled these criteria was included in analysis, 
including the control condition/s in any study exploring the in- 
teraction between active tDCS and a secondary drug or device. 
Included studies were grouped according to stimulation density, 
electrode montage, and stimulation duration (Table 1: see Sup- 
plementary material for complete MEP amplitude study selection 
flow-chart). 

As there is a growing body of research which suggests cognitive 
and/or behavioral activities undertaken during tDCS can influence 
the effects of stimulation (for discussion: Horvath et al., 2014), it is 
worth noting all studies included in the MEP analysis asked par- 
ticipants to sit comfortably and relax during stimulation. In addi- 
tion, 5 studies have explored the effect of high density tDCS on 
MEP amplitude in muscles under voluntary isometric contraction. 
Unfortunately, as two reported SDs above 100% (failed to meet 
inclusion criteria), one utilized a shoulder reference electrode lo- 
cation, and one utilized only anodal stimulation, we were unable 
to pool and analyze data from this work (see Supplementary 
material for study selection flow-chart). Finally, two studies have 
explored the effects of Hi-Definition tDCS on MEP Amplitude 
(Caparelli-Daquer et al., 2012; Kuo et al., 2013). As both of these 
came from the same research group (and utilized equipment sig- 
nificantly different from other included studies), we did not in- 
clude these in the analysis. 

As the effects of tDCS may differ following short- and long- 
duration tDCS protocols (outlined above), analyses for this mea- 
sure were initially divided into two groups: short-duration sti- 
mulation (1-7 min), and long-duration stimulation ( > 7 min). 
Additionally, several studies included in this section utilized a 
veiy-short duration stimulation protocol (4 s). As it is uncertain 
how these protocols affect outcome measures, we decide to ana- 
lyze these studies separately (veiy-short duration: < 1 min). 

Finally, to determine if technological and methodological im- 
provements have impacted the modulation effects of tDCS since 
2000, we undertook a secondary analysis by dividing papers 
which explore the impact of low-density tDCS MEP amplitude into 
three temporal cohorts {Cohort A: 2000-2004; Cohort B: 2005- 
2009; Cohort C: 2010-Present) and comparing each. This was only 
possible for this outcome measure as there were not enough pa- 
pers to undertake this analysis for any other outcome measure. 



2.3. TMS: additional measures 

See TMS method (above) for literature search and absolute 
inclusion criteria. Any study which fulfilled these criteria was in- 
cluded in analysis, including the control condition/s in any study 
exploring the interaction between active tDCS and a secondary 
drug or device. Included studies were grouped according to TMS 
outcome measure, stimulation density, electrode montage, and 
stimulation duration. Due to a lack of comparable studies, both 
low- and high-density stimulation intensities were included in a 
singular measure for cSP (Table 10: see Supplementary material 
for complete Additional TMS Measures study selection flow-chart). 
In addition, whereas two studies have explored additional TMS 
measures in muscles under voluntary isometric contraction (Jef- 
ferson et al., 2009; Kidgell et al. 2013a), one did not include nu- 
merical data; accordingly, we were unable to pool these papers. As 
all the studies that met our inclusion criteria targeted intrinsic 
hand muscles, we did not include this as a column in the relevant 
table below. Finally, again, it is worth noting all studies included in 
this section asked participants to sit comfortably and relax during 
tDCS (no cognitive or behavioral activities were undertaken during 
stimulation). 

2.4. Event related potentials (ERPs) 

See TMS method (above) for literature search and absolute 
inclusion criteria. Any study which fulfilled these criteria was in- 
cluded in analysis, including the control condition/s in any study 
exploring the interaction between active tDCS and a secondary 
drug or device. Included studies were grouped according to ERP 
component and, when possible, stimulation density, electrode 
montage, and stimulation duration. Due to a lack of comparable 
studies, varied density and/or varied electrode montages were 
pooled for several ERP measures (Table 14: see Supplementary 
material for complete ERP study selection flow-chart). With the 
exception of one study (Accornero et al., 2007), all studies included 
in this section asked participants to sit comfortably and relax 
during tDCS (no cognitive or behavioral activities were undertaken 
during stimulation). In Accornero et al. (2007), visual evoked po- 
tentials were elicited and measured during stimulation. 

2.5. EEC power spectrum 

See TMS method (above) for literature search and absolute 
inclusion criteria. Any study which fulfilled these criteria was in- 
cluded in analysis, including the control condition/s in any study 
exploring the interaction between active tDCS and a secondary 
drug or device. Included studies were grouped according to ac- 
companying behavioral activity and, when possible, electrode 
montage and stimulation duration. Due to a lack of comparable 
studies, both low- and high-density stimulation intensities were 
included for each EEC power spectrum measure (Table 19: see 
Supplementary material for complete EEC power spectrum study 
selection flow-chart). As the majority of included studies simply 
provide a written description of any significant (or non-significant) 
differences without accompanying numerical data, we were un- 
able to complete a mathematical analysis for any included mea- 
sure. Instead, we have undertaken a qualitative comparison of 
relevant studies to ascertain if an effect has reliably been found. As 
above, with the exception of one study (Wirth et al., 2011), all 
studies included in this section asked participants to sit comfor- 
tably and relax during tDCS (no cognitive or behavioral activities 
were undertaken during stimulation). In Wirth et al. (2011), a 
semantic-interference task was undertaken during stimulation. 



Please cite this article as: Horvath, J.C, et al.. Evidence that transcranial direct current stimulation (tDCS) generates little-to-no reliable 
neurophysiologic effect beyond MEP amplitude.... Neuropsychologia (2014), htlp://dx. doi. org/10. 1016/j. neuropsychologia. 2014.1 1.021 



ARTICLE IN PRESS 



].C. Horvath et al. / Neuropsychologia I fl 



2.6. Functional magnetic resonance imaging (fMRI) 

See TMS method (above) for literature search and absolute 
inclusion criteria. Any study which fulfilled these criteria was in- 
cluded in analysis, including the control condition/s in any study 
exploring the interaction between active tDCS and a secondary 
drug or device. Included studies were grouped according to elec- 
trode location, neural region of interest (ROI) analysis reported, 
and, when possible, stimulation duration. Due to a lack of com- 
parable studies, both low- and high-density stimulation intensities 
were included for each fMRI measure (Table 21 : see Supplemen- 
tary material for complete fMRI study selection flow-chart). As the 
majority of studies do not include discernable baseline measure 
values, it would be misleading to pool much of this data. In ad- 
dition, as the included studies range from 2001 to 2013, the type of 
data reported is quite variable (as later studies utilize more 
modern methods of analysis than earlier studies). Accordingly, we 
have undertaken a qualitative comparison of relevant studies to 
ascertain if an effect has reliably been found. With the exception of 
three studies (Antal et al., 2011 ; Kwon and Jang, 2011 ; Saiote et al.. 
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2013), all studies included in this section performed tDCS outside 
the fMRI and asked participants to sit comfortably/relax during 
tDCS (no cognitive or behavioral activities were undertaken during 
stimulation). In Antal et al. (2011), tDCS was delivered inside the 
scanner whilst the participants were at rest. In Kwon and Jang 
(2011), tDCS was delivered inside the scanner whilst the partici- 
pants undertook a timed grasp/release motor activity during sti- 
mulation. Similarly, in Saiote et al. (2013) tDCS was delivered in- 
side the scanner whilst the participants undertook a visuomotor 
learning task during stimulation. 



3. Results 

3.1. TMS: MEP Amplitude 



3.1.3. Low-density: primary analysis. 

3.1.1.1. Veiy-short duration stimulation (<5s). MEPs were mea- 
sured in isolation (baseline) and immediately following short 
bursts of tDCS. Modulation was determined via normalizing MEP 
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Fig. 1. Effects of tDCS on MEP amplitude (results pooled across all studies). (A) Low-density, very-short-duration ( < 1 min) tDCS; (B) low-density, short duration (1-7 min) 
tDCS; (C) low-density, long duration ( > 7 min) tDCS; (D) high-density, long-duration tDCS. Error bars: + 1 SD. 
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Table 2 

Numerical values for low-density, very-short duration tDCS on MEP amplitude. 



Table 4 

Numerical values for low-density, long duration tDCS on MEP amplitude. 



N Mean (% of SD 
baseline) 



95% u confidence interval 



Low density (0.02857 mA/cm )/very-short duration ( < 1 min) muscle at 
rest 

TO** Cathode 76 -16.34 15.26 - 19.771 < u <- 12.909 

Anode 76 26.45 30.70 19.548 < u < 33.352 

Sham - - - - 

^Difference between either stimulation condition and sham. 
Difference between stimulation conditions. 



amplitude following stimulation to baseline levels. 5 papers 
comprising 7 studies utilized this paradigm and were pooled 
(Table 1). No studies utilized a control (sham) condition. Results 
suggest that cathodal stimulation diminishes and anodal stimu- 
lation increases MEP amplitude approximately 16% and 26%, re- 
spectively (Fig. la). An unpaired student's t-test (t(150) = 10.88, 
p < 0.001) showed this difference to be significant. (Table 2). 

3.1.1.2. Short-duration stimulation (1-7 min). MEPs were measured 
pre-stimulation and continuously for 5-20 min post-stimulation. 
We pooled data at 5 min increments up to 15 min post-stimula- 
tion. 13 papers comprising 19 studies utilized this paradigm and 
were pooled (Table 1). No studies utilized a control (sham) con- 
dition. Results suggest that cathodal stimulation diminishes and 
anodal stimulation increases MEP amplitude immediately follow- 
ing stimulation approximately 27% and 38%, respectively (Fig. lb). 
At TO, an unpaired student's t-test revealed a significant difference 
between stimulation conditions (t(321) = 18.147, p< 0.001). Sig- 
nificant differences occurred between stimulation conditions at all 
time points post-stimulation (Table 3). 

3.1.1.3. Long- duration stimulation ( > 7 min). MEPs were measured 
pre-stimulation and in blocks every 5 min post-stimulation. For 
our analysis, we pooled data at 5 min increments up to 30 min 
post-stimulation and at 60 min post-stimulation. 35 papers com- 
prising 50 studies utilized this paradigm and were pooled (Ta- 
ble 1 ). 8 studies utilized a control (sham) condition - however, due 
to varying protocols, pooled sham data could only be obtained 
immediately, 10- and 30 min post stimulation. Results suggest that 
cathodal stimulation diminishes and anodal stimulation increased 
MEP amplitude immediately following stimulation approximately 
26% and 35%, respectively (Fig. Ic). A one-way ANOVA (including 

Table 3 

Numerical values for low-density, short duration tDCS on MEP amplitude. 



Mean 
(% of 
baseline) 



SD 



95% II confidence interval 



Low density (0.02857 niA/cm^)/short duration (1-7 min) muscle at rest 



TO** 


Cathode 


173 


-27.29 


19.93 


-30.260<u< - 24.320 




Anode 


150 


3754 


41.84 


30.844 <u<44.236 




Sham 










15** 


Cathode 


152 


-14.20 


22.14 


- 17720 < u < - 10.680 




Anode 


141 


24.54 


3737 


18.372 <u < 30.708 




Sham 










110** 


Cathode 


173 


-5.84 


21.62 


- 9.062 <u < -2.618 




Anode 


158 


11.92 


31.01 


7.085 <u < 16.755 




Sham 










115** 


Cathode 


103 


-6.18 


20.57 


-10.153<u< - 2.207 




Anode 


102 


710 


19.36 


3.343 < u < 10.857 




Sham 















N 


Mean (% of 
baseline) 


SD 


95% u confidence interval 


Low density (0.02857 mA/cm^)/long duration ( > 7 min) muscle at rest 


TO* 


Cathode 


400 


-25.05 


28.82 


-27.883 <u< -22.217 




Anode 


412 


35.08 


46.92 


30.549 < 11 < 39.611 




Sham 


64 


1.81 


20.82 


-3.291 <u < 6.911 


T5** 


Cathode 


287 


-29.44 


21.98 


-31.983 <u < -26.897 




Anode 


320 


32.01 


35.42 


28.129 <u < 35.891 




Sham 


- 


- 


- 


- 


TIO* 


Cathode 


284 


-29.26 


25.81 


-32.262 <u< -26.258 




Anode 


341 


32.88 


3778 


28.870<u<36.890 




Sham 


32 


-3.88 


30.82 


- 14.559 <ti< 6.799 


T15** 


Cathode 


227 


-26.79 


1785 


-29.112 <u< - 24.468 




Anode 


290 


30.99 


35.86 


26.863 <u<35.117 




Sham 










T20** 


Cathode 


262 


-26.24 


19.51 


-28.613<u< -23.867 




Anode 


290 


28.40 


32.76 


24.629 <u< 32.171 




Sham 










T25* 


Cathode 


227 


-26.75 


1738 


-29.011 <u< -24.489 




Anode 


290 


25.88 


33.13 


22.067 <u< 29.693 




Sham 










T30* 


Cathode 


289 


-24.14 


26.10 


-27.149 <u< - 21.131 




Anode 


299 


22.87 


32.10 


19.231 < u < 26.509 




Sham 


35 


1.6 


51.40 


- 15.429 <ii< 18.629 


T60* 


Cathode 


220 


-19.68 


19.71 


-22.285 <u< -17.075 




Anode 


274 


13.15 


23.48 


10.370 < ti< 15.930 




Sham 











^Difference between either stimulation condition and sham. 
Difference between stimulation conditions. 



' Difference between either stimulation condition and sham. 
Difference between stimulation conditions. 



sham) at TO revealed a significant difference between conditions (F 
(2, 873)=253.78, p< 0.001). Post-hoc tests revealed a significant 
difference between cathode and anode (f(810)=21.93, p< 0.001), 
between cathode and sham (t( 462)= 7.16, p < 0.001), and between 
anode and sham (t(474)=5.58, p < 0.001 ).Significant differences 
occurred between stimulation conditions at all time points post- 
stimulation and between stimulation conditions and sham at TIO 
and T30 (Table 3 and 4). 

3.1.2. High-density 

3.1.2.1. Long-duration stimulation (>7min) - tested muscle at 
rest MEPs were measured pre-stimulation and in blocks eveiy 5 min 
post-stimulation. For our analysis, we pooled data at 5 min increments 
up to 30 min post-stimulation and at 60 min post-stimulation. 13 
papers comprising 16 studies utilized this paradigm and were pooled 
(Table 1 ). 5 studies utilized a control (sham) condition - however, due 
to varying protocols, pooled sham data could only be obtained im- 
mediately post stimulation, at T30, and at T60. Results suggest that 
cathodal stimulation diminishes and anodal stimulation increased 
MEP amplitude immediately following stimulation approximately 22% 
and 30%, respectively (Fig. Id). A one-way ANOVA (including sham) at 
TO revealed a significant difference between conditions (f(2, 324)= 
33.228, p < 0.001). Post-hoc tests revealed a significant difference be- 
tween cathode and anode (t(272)= 7.621, p < 0.001 ), between cathode 
and sham (£(183)=2.438, p=0.016), and between anode and sham (f 
(193) =3.868, p< 0.001).. Significant differences occurred between 
stimulation conditions at all measurable time points (Table 5). Inter- 
estingly, corrected t-tests at T30 and T60 reveal that, although anode 
differs from sham at both these time points (T30: t( 174) =4.797, 
p< 0.001; T60: t(125)= 2.590, p= 0.011), cathode does not (T30: f 
( 107) = 1.344, p = 0.181; T60: t(123)= 1.102, p= -0.273). 

3.1.3. Low-density: secondary analysis (temporal cohorts). 

3.1.3.1. Very-short duration stimulation (<1 min). As the last time 
this measure was reported in the literature was 2007, only Cohorts 
A and B could be compared. An unpaired student's t-test revealed 
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Table 5 

Numerical values for high-density, long duration tDCS on MEP amplitude. 







JV 


Mean (% of 
baseline) 


SD 


95% u confidence interval 


High density ( > 0.02857 mA/cm^)/long duration ( > 7 min) - muscle at rest 


TO* 


Cathode 


132 


-11.66 


33.87 


-17.491 <u< - 5.828 




Anode 


142 


2713 


48.51 


19.083 <u < 35.178 




Sham 


53 


0.49 


20.36 


-5.121 <u < 6.101 


T5** 


Cathode 


32 


-29.25 


19.83 


-36.121 <u< - 22.379 




Anode 


54 


20.04 


43.97 


12.885 <u<32.195 




Sham 


- 


- 


- 


- 


TIO** 


Cathode 


40 


-20.55 


25.38 


-28.667 <u< -12.433 




Anode 


71 


23.72 


34.52 


15.549 < IK 31.891 




Sham 


- 


- 


- 


- 


T15** 


Cathode 


48 


-23.64 


24.62 


-30.789<u< -16.491 




Anode 


70 


35.54 


49.93 


23.635 <u<47.445 




Sham 










T20** 


Cathode 


39 


-35.60 


29.57 


-44.881 <u< - 26.319 




Anode 


54 


34.37 


42.17 


23.122 <u<45.618 




Sham 










T25 


Cathode 












Anode 


36 


2733 


32.55 


32.5531 




Sham 










T30* 


Cathode 


77 


-6.20 


45.00 


-16.414<u< 4.014 




Anode 


99 


26.27 


44.20 


17.454 < u < 35.086 




Sham 


32 


4.97 


20.63 


-2.468 <u< 12.408 


T60* 


Cathode 


93 


-1.65 


33.64 


-8.578<u<5.278 




Anode 


95 


29.17 


47.65 


19.463 <u< 38.877 




Sham 


32 


5.81 


31.13 


-5.414<u< 17034 



' Difference between either stimulation condition and sham. 
' Difference between stimulation conditions. 



no significant difference in reported catiiodal effects (t(74)= 1.490, 
p= 0.140). However, an unpaired student's t-test revealed a sig- 
nificant difference in reported anodal effects (t(74) = 5.265, 
p < 0.001 ) sucfi tliat tiie effects are significantly smaller in more 
recent studies (Table 6; Fig. 2b). 



3.1.3.2. Short-duration stimulation (1-7 min). One-way ANOVAs 
were conducted at each time point pooled post-stimulation (Ta- 
ble 7). Of note, recent studies (Cohort C) showed significantly 
smaller effects immediately following anodal stimulation com- 
pared to earlier studies (Fig. 2c,d). 



3.1.3.3. Long-duration stimulation (> 7 min). One-way ANOVAs 
were conducted at each time point pooled post-stimulation (Ta- 
ble 8). Of note, recent studies (Cohort C) showed significantly 
smaller effects immediately following cathodal stimulation com- 
pared to earlier studies {Cohort A). In addition. Cohort C showed 
significantly smaller effects at every time point post-anodal sti- 
mulation compared to Cohorts A and B (Fig. 2e and f). 



3.2. TMS: additional measures 

3.2.1. Resting and active motor threshold (rMT and aMT) intensity 
rMTs were measured prior to and immediately following sti- 
mulation only. 6 studies have explored the impact of tDCS on this 
measure and were pooled for rMT analysis (Table 9). 1 study uti- 
lized a control (sham) condition. A one-way ANOVA revealed no 
significant difference between conditions (F(2, 115)=0.056, 
p=0.946: Table 11; Fig. 3a). Table 10). An additional 8 studies 
which have explored the impact on this measure reported no ef- 
fect - although none supplied numerical data (see Table 10 notes). 

aMTs were measured prior to and immediately following sti- 
mulation only. 5 studies have explored the impact of tDCS on this 
measure and were pooled for aMT analysis (Table 9). A one-way 
ANOVA revealed no significant difference between conditions (f(2, 
117)=0.585, p=0.559: Table 10; Fig. 3b). 

4 studies from 3 papers have explored the effects of veiy-short 
duration ( < 5 sec) and short-duration (1-7 min) stimulation on 
rMT (Cengiz et al., 2013; Mordillo-Mateos et al., 2012; Nitsche 
et al., 2005). Although each reported no effect of tDCS on rMT, due 
to varying protocols, no data could be pooled. 

3.2.2. Cortical silent period (cSP) duration 

cSPs were measured prior to and immediately following sti- 
mulation only. 7 studies have explored the impact of tDCS on this 
measure. Each reported no effect of tDCS on cSP. However, only 
3 of these studies reported numerical data and could be pooled 
(Table 9). One study utilized a control (sham) condition. A one-way 
ANOVA revealed no significant difference between conditions (f(2, 
53)=0.751, p=0.477: Table 11; Fig. 3c). 

3.2.3. Short intei-val cortical inhibition (SICI) and intracortical facil- 
itation (ICF) duration 

For this analysis, we combined all ISIs thought to lead to SlCl (2, 3, 4, 
and 5 ms) into a single measure and all ISIs thought to lead to ICF (7, 9, 
10, 12,and 15ms) into a single measure. So as not to compound N, 
weighted averages and pooled standard deviations were determined for 
any paper that explored more than one of these ISIs. In addition, a post- 
stimulation MEP amplitude comparison was made between groups. 

SICI and ICF were typically measured prior to and immediately 
following stimulation only. 8 studies have explored the impact of 
tDCS on these measures and were pooled for analysis (Table 9). 
Two studies utilized a control (sham) condition. With regards to 
SICI, student's f-tests (Cathode: t(172)=0.076, p=0.940; Anode: t 
(82) = 1.616, p = 0.110; Sham: t(30) = 0.426, p = 0.673; Fig. 3d) did 
not reveal a significant difference between pre- and post-scores for 
any condition (Table 12). A one-way ANOVA between post-sti- 
mulation scores did reveal a significant difference (f(2, 122)= 
9.717, p < 0.001 ). Post-hoc bonferroni correct t-tests revealed no 
difference between cathode and sham (t(81)= 1.477, p=0.144) or 
anode and sham (f(56) = 1.420, p = 0.1611), but a significant dif- 
ference between anode and cathode (t(107)=4.931, p < 0.001). 



Table 6 

Temporal analysis numerical values for low-density, very-short duration tDCS on MEP amplitude. 



67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
91 
92 
93 
94 
95 
96 
97 
98 
99 
100 
101 
102 
103 
104 
105 
106 
107 
108 
109 
110 
111 
112 
113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 
131 
132 



N 



Mean (% of baseline) 



SD 



95% u confidence interval 



Low density (0.02857 mA/cm )/very-short duration ( < 1 min) - muscle at rest 

Cathode TO 2000-2004 40 -18.85 

2005-2009 36 -13.67 

2010-Present 

Anode TO* 2000-2004 40 41.55 

2005-2009 36 9.67 
2010-Present 



6.95 
20.75 



32.40 
17.30 



-21.004<u< -16.696 
-20.448 <u < -6.892 



31.509 <u < 51.591 
4.019 <u < 15.321 



' Difference between temporal cohorts. 
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Fig. 2. Effects of tDCS on MEP amplitude (results pooled according to temporal cohorts). (A and B) Low-density, very-short-duration ( < 1 min) cathodal and anodal tDCS. 
respectively; (C and D) low-density, short duration (1-7 min) cathodal and anodal tDCS, respectively; (E and F) low-density, long duration ( > 7 min) cathodal and anodal 
tDCS. respectively. Error bars: + 1 SD. 
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Table 7 

Temporal analysis numerical values for low-density, short duration tDCS on MEP amplitude. 







N 


Mean (% of baseline) 


SD 


95% u confidence interval 


Low density (0.02857 mA/cm^)/short duration (1-7 min) 


— muscle at rest 










Cathode TO 


2000-2004 


90 


-31.34 




16.57 


-34.763 <u< -27917 




2005-2009 


48 


-25.75 




11.63 


- 29.040 < u < -22.460 




2010-Present 


26 


-16.54 




29.83 


-28.591 <u< -4.489 


T5* 


2000-2004 


69 


-7.71 




2236 


-13.033 < u < -2.387 




2005-2009 


48 


-23.50 




12.17 


-27093 <u< -19.907 




2010-Present 


26 


-8.54 




31.07 


-21.092 <u<4.012 


no* 


2000-2004 


90 


-1.72 




22.80 


-6.431 <u<2.991 




2005-2009 


48 


- 16.00 




13.44 


-19.802 < u < -12.198 




2010-Present 


26 


-2.62 




26.03 


- 13.136 < u < 7.896 


T15 


2000-2004 


30 


0 00 




25.25 


^ - 9.036 <ti< 9.036 




2005-2009 


48 


-10.5 




14.53 


-14.611 <u< -6.389 




2010-Present 


16 


-6 




28 71 




Anode TO* 


2000-2004 


64 


44.67 




43.73 


33.956<u< 55.384 




2005-2009 


48 


36.00 




22.75 


29.564<u< 42.436 




2010-Present 


46 


25.61 




26.65 


17.909 <u<33.311 


T5 


2000-2004 


57 


18.74 




^' 40.88 


8.127 <u<29.353 




2005-2009 


48 


22.00 




16.18 


17423 <u< 26.577 




2010-Present 


36 


19.11 




30.90 


9.016 < u < 29.204 


TIO* 


2000-2004 


64 


5.09 




16.89 


0.952 <u< 9.228 




2005-2009 


48 


15.00 




16.11 


10.442 <u< 19.558 




2010-Present 


38 


15.95 


V 


39.81 


3.292 < u < 28.608 


T15 


2000-2004 


18 


1.33 


16.04 


-6.647<u<9.307 




2005-2009 


48 


8.00 


15.21 


3.369 < u < 12.303 




2010-Present 


28 


9.86 




27.46 


-0.830<u<20.550 



67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
91 
92 
93 
94 
95 
96 
97 
98 
99 
100 
101 
102 
103 
104 
105 
106 
107 
108 
109 
110 
111 
112 
113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 
131 
132 



' Difference between temporal cohorts. 



With regards to ICF, student's t-tests (Cathode: t(132)=0.549, 
p=0.584; Anode: t(82)=0.799, p=0.427; Sham: t(24)= 1.235, 
p=0.229; Fig. 3e) did not reveal a significant difference between 
pre- and post-scores for any condition (Table 12). A one-way AN- 
OVA between post-stimulation scores only did not reveal a sig- 
nificant difference between groups (F(2,122)= 1.700, p = 0.187). 

3 studies have examined the effect of very-short duration 
( < 5 s) and short-duration (1-7 min) stimulation on SICI and ICF 
((Fricke et al., 2011; Munneke et al., 2011; Nitsche et al., 2005). 
However, as one of these studies reported SDs far above 100% and 
the remaining two came from the same research group, we did not 
pool this data for analysis. 

3.3. Event related potentials (ERPs) 

3.3.1. PlOO visual evoked potential (VHP) amplitude: low and high 
contrast stimuli 

PlOO amplitudes were typically measured pre-stimulation and 
in blocks every 10 min post-stimulation (up to 30 min). 3 Studies 
met criteria for inclusion in the PlOO analysis (Table 13). No stu- 
dies utilized a sham control condition. With regards to low-con- 
trast stimuli, student's f-tests revealed no significant difference 
between stimulation conditions at any time point (TO: f(50) = 
1.472, p = 0.147; TIO: t(50) = 1.768, p = 0.084; T20: t(50) = 0.779, 
p=0.440; T30: t(50)=0.492, p=0.625; Fig. 4a; Table 14). With 
regards to high contrast stimulation, student's t-tests revealed no 
significant difference between stimulation conditions at any time 
point (TO: t(61)=0.931,p=0.356; TIO: t(50) = 0.816,p = 0.419; T20: 
t(50) = 0.080, p=0.937; T30: t(50)=0.837, p = 0.406; Fig. 4b; Ta- 
ble 14) In addition, all included studies reported no effect of tDCS 
on the latency of the PlOO component; although no discernable 
numerical data was available for pooling. 

3.3.2. N20 somatosensory evoked potential (SEP): median nerve 
electrical stimuh 

N20 amplitudes were measured pre-, immediately, 15 min, and 
20 min post-stimulation. Due to constraints in the sham condition, 



time points T15 and T20 were combined into a single measure for 
each group (whilst maintaining N). 3 Studies met criteria for in- 
clusion in this analysis (Table 13). One study utilized a control 
(sham) condition. A student's t-tests at TO revealed no significant 
difference between stimulation conditions (TO: F(2, 67)= 1.895, 
p=0.158; T20: f(2)=1.220, p=0.302). This non-significance con- 
tinued at each time point post-stimulation (Table 15; Fig. 4c). In 
addition, all included studies reported no effect of tDCS on the N30 
SEP component; though, no discernable numerical data was 
available for pooling. 

3.3.3. N2 and P2 laser evoked potential (LEP): heat pain stimuli 
N2 and P2 amplitudes were measured prior to and immediately 

following stimulation only. 3 Studies met criteria for inclusion in 
this analysis (Table 13). One-way ANOVAs revealed no significant 
difference between groups for either outcome measure (N2: f(2, 
76) = 0.625, p = 0.538; P2: f(2, 76) = 0.997, p = 0.374; Table 16; 
Fig. 4d and e). In addition, all included studies noted no effect of 
tDCS on the amplitude of the Nl LEP component or any compo- 
nential latency; although not enough data was available for 
pooling. Finally, (Terney et al., 2008) and (Csifcsak et al., 2009) 
further noted no effect of tDCS on LEP threshold; though, again, no 
discernable numerical data was available for pooling. 

3.3.4. Mu/Alpha event related desynchronization (ERD): motor 
imagery 

Mu/Alpha ERD percentage was measured prior to and im- 
mediately following stimulation only. 3 studies met criteria for 
inclusion in this analysis (Table 13). A one-way ANOVA revealed no 
significant difference between groups (F(2, 94) = 0.420, p=0.658; 
Table 17; Fig. 4f). 

3.4. EEC power specti-um results 

3.4.1. Oscillation frequency power: at rest 

4 studies have explored the effect of tDCS on resting state EEC 
oscillatory power using comparable protocols (Table 18). With 
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Table 8 

Temporal analysis numerical values for low-density, long duration tDCS on MEP amplitude. 







N 




Mean (% of baseline) 


SD 


95% u confidence interval 


Low density (0.02857 niA/cm^)/long duration ( > 7 min) - 


muscle at rest 










Cathode TO* 


2000-2004 


108 




-33.60 


21.37 


-37.630 < u < -29.570 




2005-2009 


128 




-22.25 


22.12 


-26.082 < u < -18.418 




2010-Present 


164 




22.13 


35.90 


- 27.665 < u < - 16.595 


T5* 


2000-2004 


72 




-36.64 


21.54 


-41.615 <u< -31.665 




2005-2009 


132 




-25.87 


19.92 


-29.268 < u < -22.472 




2010-Present 


71 




-29.35 


24.76 


-35.109 < u < -23.591 


TIG 


2000-2004 


72 




-32.75 


18.25 


-36.996 < u < -28.534 




2005-2009 


132 




-27.89 


19.41 


-31.201 < u < -24.579 




2010-Present 


62 




-24.35 


37.98 


-33.804 < u < - 14.896 


T15* 


2000-2004 


72 




-32.19 


15.96 


-35.877 < u < -28.503 




2005-2009 


122 




-24.14 


18.02 


-27.338 < u < -20.942 




2010-Present 


33 




-24.82 


19.00 


-31.303 < u < -18.337 


T20* 


2000-2004 


72 




-31.92 


16.44 


-35.717<u< -28.123 




2005-2009 


122 




-24.90 


17.38 


-27.984 < u < -21.816 




2010-Present 


68 




-22.60 


26.50 


-29.014<u< -16.186 


T25 


2000-2004 


72 




-29.75 


19.63 


-34.284<u< -25.216 




2005-2009 


122 




-26.54 


15.36 


-29.266<u< -23.814 




2010-Present 


33 




-21.00 


18.23 


-27220 < u < -14.780 


T30 


2000-2004 


72 




-27.39 


16.96 


-31.308 < u < -23.472 




2005-2009 


122 




-25.48 


19.42 


-28.926<u< -22.034 




2010-Present 


95 




-20.18 


36.83 


-27586 < u < -12.774 


T60 


2000-2004 


102 




-18.82 


17.63 


-22.241 <u< -15.399 




2005-2009 


125 




-21.15 


19.63 


-24.591 < u < -17709 




2010-Present 


42 




- 23.07 


20.02 


-29.125 <u < -17015 


Anode TO* 


2000-2004 


123 




50.87 


44.78 


42.956 < u < 58.784 




2005-2009 


128 




36.38 


37.20 


29.935 < u < 42.825 




2010-Present 


176 




18.65 


50.94 


11.124 <ii< 26.176 


T5* 


2000-2004 


87 




51.09 


41.07 


42.260 < u < 59.720 




2005-2009 


132 




32.58 


28.64 


27694 < u < 37466 




2010-Present 


96 




9.019 


31.05 


2.808 < u< 15.230 


TIG* 


2000-2004 


87 


■I 


47.33 


38.53 


39.234 < u < 55.426 




2005-2009 


132 




30.95 


27.07 


26.332 <u <35.568 




2010-Present 


137 




19.70 


44.14 


12.309 < u < 27091 


T15* 


2000-2004 


87 




48.31 


37.61 


40.407 < u < 56.213 




2005-2009 


122 




32.88 


32.91 


27040 < u < 38.720 




2010-Present 


96 




4.29 


28.94 


-1.499 <u< 10.079 


T20* 


2000-2004 


87 




34.01 


30.25 


27653 < u < 40.367 




2005-2009 


122 




35.32 


34.35 


29.225 < u < 41.415 




2010-Present 


96 




6.47 


29.70 


0.529 < ti< 12.411 


T25* 


2000-2004 


87 




38.39 


32.41 


31.580 <u<45.200 




2005-2009 


122 




32.57 


29.64 


27310 < u < 37830 




2010-Present 


96 




-1.30 


27.36 


-6.733 <u<4.173 


130* 


2000-2004 


87 




32.34 


28.66 


26.318 <u<38.362 




2005-2009 


122 




25.78 


29.53 


20.540 < u < 31.020 




2010-Present 


105 




6.08 


35.17 


-0.647<u< 12.807 


T60* 


2000-2004 


77 




16.55 


28.05 


10.285 <u<22.815 




2005-2009 


116 




16.79 


21.09 


12.952 <u<20.628 




2010-Present 


96 




1.48 


22.25 


-2.971 <u< 5.931 



' Difference between temporal cohorts. 



regards to cathodal stimulation, (Ardolino et al., 2005) reported a 
significant increase in the power of delta oscillations (compared to 
sham), (Pellicciari et al., 2013) reported an increase in the power of 
alpha and theta oscillations following cathodal stimulation (al- 
though this increase was comparable to increases following anodal 
stimulation and no sham condition was used), (Notturno et al., 
2013) reported a significant decrease in power of alpha oscillations 
(compared to sham), and (Matsumoto et al., 2010) reported no 
effect on any oscillatory frequency (compared to sham). Each pa- 
per reported no significant effect of cathodal stimulation on any 
other oscillatory frequency (Table 19). 

With regards to anodal stimulation, (Pellicciari et al., 2013) 
reported an increase in the power of alpha and theta frequencies 
(although this increase was comparable to increases following 
cathodal stimulation and no sham condition was used), whilst 
(Matsumoto et al., 2010) and (Notturno et al., 2013) both reported 
no effect of anodal stimulation on any oscillatory frequency 
(compared to sham: Table 19). 



3.4.2. Oscillation frequency power (and P300 ERP): during N-Back 
working memory task 

2 Studies have explored the effect of tDCS on EEC oscillatory power 
and P300 ERP characteristics during an N-Back task using comparable 
protocols (Table 18); unfortunately, as one of these studies only ex- 
plored anodal stimulation, we will be unable to discuss the effects of 
cathodal stimulation. (Keeser et al., 2011) reported a reduction of ab- 
solute Delta power across frontal electi'odes during a 2-back task (but 
not during a 0- or 1-back task: compared to sham). These researchers 
did not find any effect of tDCS on the theta, alpha, beta, or gamma 
bands (for 10 min following stimulation) during any N-Back condition. 
In addition, this group reported no significant change in P300 ampli- 
tude or latency at Pz (although they reported some attenuation at 
other electrodes, Pz was the only electrode explored in the second 
included study). Zaehle et al. (2011) reported no significant differences 
in frequency power following anodal stimulation (compared to sham). 
In addition, these researchers reported no significant difference in 
P300 amplitude or latency at Pz (Table 19). 
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Table 9 

Studies exploring tlie effects of long duration ( > 7 min) tDCS on TMS measured resting motor tfireshold levels, active motor threshold levels, cortical silent period, short 
interval cortical inhibition, and intracortical facilitation. 



Paper/# of studies 


N 


tDCS duration 


tDCS polarities 


DB/NN 


Parameters 


#Averaged 


Resting motor threshold 














Lang et al. (2004) 


5 


10 min 


A, C 


1 


5 of 10 MEPs between 50-100 uV 




Quartarone et al. (2005) 


8 


10 min 


A, C, S 


-/- 


5 of 10 MEPs above 50 uV 




Nitsche et al. (2005) 


12 


9 min C/11 min A 


A, C 


-/- 


3 of 6 MEPs above 50 uV 




Scelzo et al. (2011) 


12 


13 min 


A 




5 of 10 MEPs above 50 uV 




Hasan et al. (2012) 


18 


9 min 


C 




5 of 10 MEPs above 50 uV 




Di Lazzaro et al. (2012) 


30 


20 min 


C 




5 of 10 MEPs above 50 uV 




Active motor threshold 














Lang et al. (2004) 


5 


10 min 


A, C 




5 of 10 MEPs between 200-300 uV (contraction) 




Quartarone et al. (2005) 


8 


10 min 


A, C, S 




5 of 10 MEPs above 200 uV (contraction) 




Nitsche et al. (2005) 


12 


9 min C/11 min A 


A, C 




3 of 6 MEPs greater than background activity (15% 
contraction) 




Thirugnanas-ambandam et al. 


16 


20 min 


A, C 




5 of 10 MEPs above 100 uV (peak-to-peak: 15% 




(2011) 










contraction) 




Di Lazzaro et al. (2012) 


30 


20 min 


C 




5 of 10 MEPs above 200 uV (20% contraction) 




Cortical silent period (cSP) 














Hasan et al. (2012) 


18 


9 min 0.258 mA/cnr' 


C 


-1- 


1 mV MEP w/25-30% contraction 


10? 


Suzuki et al. (2012) 


6 


10 minO.258 wA/cm^ 


A, C, S 


-1- 


110% rMT w/10-20% contraction 


WRel cSP 


Tremblay et al. (2013) 


10 


20 min 0.6 mA/crn^ 


A, C 


-/Y 


1 mV MEP w/25-30% contraction 


WAbs cSP 


Short interval cortical inhibition (SICI) and intracortical facilitation (ICF) 








Siebner et al.O (2004) 


8 


10 min 


A, C 


-1- 


Test Pulse: 1 mV MEP Cond. Pulse: 80% aMT 


15 


Nitsche et al. (2004) 


6 


11 min 


A 


-1- 


Test Pulse: 1 mV MEP Cond. Pulse: 70% aMT 


15 


Quartarone et al. (2005) 


8 


10 min 


A, C, S 


-1- 


Test Pulse: 0.8-1 mV MEPCond. Pulse: 80% aMT 


15 


Nitsche et al. (2005) 


12 


9 min C/13 min A 


A, C 


-1- 


Test Pulse: 1 mV MEP Cond. Pulse: 70% aMT 


15 


Boros et al. (2008) 


8 


13 min 


A 


-1- 


Test Pulse: 1 mV MEP Cond. Pulse: 70% aMT 


12 


Di Lazzaro et al. (2012) 


30 


20 min 


C 


-1- 


Test Pulse: 1 mV MEP Cond. Pulse: aMT - 5% MSO 


5 


Batsikadze et al. (2013/x2) 


9 (xl)8 
(xl) 


20 min (xl ) 


C (xl)S (xl) 


-1- 


Test Pulse: 1 mV MEP Cond. Pulse: 70% aMT 


15 



General notes: DB= double-blind experimental design with regards to tDCS condition (not any included pharmacological arm); NN^neuronavigation system used to guide 
TMS; r(x)=time (in minutes) post stimulation that measurements were reported and pooled for this analysis. A=anode; C=cathode; S = sham; Rel cSP=relative cortical 
silent period; Abs cSP= absolute cortical silent period. 

rMT and aMT Notes: Alonzo et al. (2012), Batsikadze et al. (2013), Cambieri et al. (2012), Hasan et al. (2013), Mordillo-Mateos et al. (2012), Schmidt et al. (2013), Simis et al. 
(2013), and Thirugnanasambandam et al. (2011) were excluded from analysis as they did not supply numerical data. 

cSP notes: Thirugnanasambandam et al. (2011), Batsikadze et al. (2013), and Hasan et al. (2013) were excluded from analysis as they did not supply discernable numerical 
data. However, each of these papers reported no significant change in cSP due to stimulation. 

SIG and ICF notes: two studies from Munneke et al. (2011 ) were excluded from analysis for reporting SDs above 100%. Thirugnanasambandam et al. (2011 ) was excluded from 
analysis for reporting only percent change (rather than percentage of unconditioned response). Simis et al. (2013) was excluded from analysis as they did not supply 
numerical data (although this paper reported no significant effect of stimulation on either SlCl or ICF). Jefferson et al. (2009) was excluded from analysis as they did not 
supply numerical data (although this paper reported no significant effect of stimulation on either SICI or ICF). 



3.4.3. Oscillation frequency characteristics: during picture tasks 

2 studies have explored the effect of tDCS on EEC oscillatory 
power during picture naming tasl<s using comparable protocols 
(Table 18); unfortunately, as both only explored anodal stimula- 
tion, we will be unable to discuss the effects of cathodal stimula- 
tion. (Wirth et al., 2011) reported no significant changes in power 
across any frequency following anodal stimulation during a neu- 
tral valence picture naming task (compared to sham). Similarly, 
(Maeoka et al., 2012) reported no significant changes in power 
across any frequency following anodal stimulation during either 
neutral or positive valence picture viewing tasks (compared to 

Table 10 

Numerical values for low-density, long duration tDCS on rMT and aMT values. 

iV Mean (% of baseline) SD 95% u confidence interval 



Resting motor threshold (rMT) 

TO Cathode 73 1.15 

Anode 37 0.22 

Sham 8 0 
Active motor threshold (aMT) 
iV Mean (% of baseline) 

TO Cathode 71 0.99 

Anode 41 - 1,76 

Sham 8 0 



16.94 - 2.736<u<5.036 

12.41 -3.779 <u<4.219 

14 -11.681 <u< 11.681 

SD 95% u confidence interval 

14.26 -2.327 <u<4.307 

10.24 -4.894<u< 1.374 

13 -10.847<u< 10.847 



sham). However, these researchers did report a significant reduc- 
tion in beta and increase in alpha band power (compared to sham) 
during an unpleasant valence picture viewing task (Table 19). 



3.5. fMRI results 

3.5.3. Short-duration stimulation (1-7 min) 

2 studies have explored the effect of very-short-duration and/ 
or short-duration stimulation over Ml on voxel activation/spread 
under the active electrode during an active motor task (Table 20). 
Although each utilized a slightly different protocol, neither re- 
ported a significant change in either measure following either 
anodal or cathodal stimulation (Table 21 ). 



Table 11 

Numerical values for long duration tDCS on cSP duration. 

IV Mean (% of baseline) SD 95% u confidence interval 



Cortical silent period (cSP) 

TO Cathode 34 3 
Anode 16 -6.75 
Sham 6 0 



27.49 
19.02 
35 



-6.240<u< 12.240 
-16.878 <u<3.378 
-36.722 <u<36.722 



"Difference between either stimulation condition and sham. 
''Difference between stimulation conditions. 



'Difference between either stimulation condition and sham. 
''Difference between stimulation conditions. 



67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
91 
92 
93 
94 
95 
96 
97 
98 
99 
100 
101 
102 
103 
104 
105 
106 
107 
108 
109 
110 
111 
112 
113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 
131 
132 
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Fig. 3. Effects of tDCS on various TMS outcome measures. (A and B) Low-density, long duration ( > 7 min) tDCS exploring resting and active motor thresholds, respectively; 
(C) long-duration tDCS exploring cortical silent period duration; (d and E) low-density, long duration tDCS exploring short interval cortical inhibition and intracortical 
facilitation, respectively. Error bars: + 1 SD. 



3.5.1.1. Long-duration stimulation (>7min). 2 studies have ex- 
plored the effect of long-duration catliodal stimulation over Ml on 
voxel activation/spread under the active electrode during an active 
motor task (Table 20). Neither reported any significant change in 



Table 12 

Numerical values for low-density, long duration tDCS on SlCl and ICF. 



N Mean (% SD 
of base- 
line) 



95% u confidence interval 



Short interval cortical inhibition (SICI) and intracortical facilitation (ICF) 



sicr2, 3, 


Cathode 


Pre 


67 


-40.69 


44.05 


-51.238 <ti< -30.142 


4, and 




Post 


67 


-46.78 


33.45 


-54.790<u< -38.77 


5 ms 


Anode 


Pre 


42 


-27.61 


49.59 


-42.608 <ti< -12.612 


ISI 




Post 


42 


-11.67 


40.32 


-23.864<u<0.524 




Sham 


Pre 


16 


-31.88 


30.65 


-48.201 <u< -15.559 






Post 


16 


-27.38 


29.13 


-42.892 <u< -11.868 


ICF 7, 9, 


Cathode 


Pre 


67 


20.93 


33.96 


12.798 <u< 29.062 


10, 12, 




Post 


67 


23.99 


30.44 


16.701 < ii< 31.279 


and 


Anode 


Pre 


42 


21.93 


40.24 


9.760 <u<34.100 


15 ms 




Post 


42 


29.12 


42.26 


16.339 <u<41.901 


IS! 


Sham 


Pre 


13 


19.69 


32.90 


-0.202 <u<39.582 






Post 


13 


44.54 


64.67 


5.439 < u < 83.641 



'Difference between pre- and post-stimulation. 

''Post-stimulation score difference between either stimulation condition and sham. 
^Post-stimulation score difference between stimulation conditions. 



either measure. 3 studies have explored these measures following 
anodal stimulation. One study reported a significant increase in 
activation and spread whilst the other two reported no significant 
change in either measure following stimulation (Table 21 ). 



3.5.2. ROI: SMA activation/spread: during active motor taslc. 
3.5.2.1. Short-duration stimulation (1-7 min). 3 Studies have ex- 
plored the effect of immediate and/or short-duration anodal sti- 
mulation over Ml on voxel activation/spread at SMA during an 
active motor task (Table 20). One study reported a significant in- 
crease, one study reported a significant decrease, and one study 
reported no significant change in both measures (Table 21). 2 stu- 
dies have explored these measures following cathodal. One study 
reported a significant decrease whilst the other reported no sig- 
nificant change in both measures (Table 21). 



3.5.2.2. Long-duration stimulation ( > 7 min ). 3 Studies have ex- 
plored the effect of long-duration anodal stimulation over Ml on voxel 
activation/spread at SMA during an active motor task (Table 20). Two 
studies reported a significant increase whilst one reported an insig- 
nificant decrease in both measures (Table 21 ). 2 studies have explored 
this measure following cathodal stimulation. Neither reported any 
significant change on either measure (Table 21). 
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Table 13 

Studies exploring tlie effects of long duration ( > 7 min) tDCS on varied event related potential (ERP) measures. 



Paper/* of studies 


N 


tDCS duration (min) 


tDCS polarities 


DB Recording electrode 


# Stimuli averaged 


Hx) 


PlOO VEP amplitude: low and high contrast stimuli 










Antal et al. (2004) 


16 


10 (xl) 


A, C 


Oz 


50 


T0-T30 


Accornero et al. (2007) 


10 


10 


A, C 


Oz 


60 


T0-T30 


Vigano et al. (2013) 


11 


15 


A 


— Oz(High Contrast only) 


600 


TO 


N20 SEP amplitude 














Dieckhofer et al. (2006) 


10 


9 


A, C 


CP3 


? 


TO, T20 


Kirimoto et al. (2011/x2) 


10 (x2) 


15 (x2) 


A, C, S (xl )A, C (xl ) 


CPS 


300 


TO, T15-20 


N2 and P2 LEP amplitude 














Antal et al. (2008) 


10 


15 


A, C, S 


Cz 


40 


TO 


Terney et al. (2008) 


12 A,C, 7 S 


15 


C, S 


Cz 


40 


TO 


Csifcsak et al. (2009) 


10 


10 


A, C, S 


Cz 


40 


TO 


Mu/alpha ERD: motor imagery 












Matsumoto et al. (2010) 


6 


10 


A, C, S 


C3 and surrounds 


n/a 


TO 


Wei et al. (2013) 


8 A. 8 S 


15 


A, S 


C3 and surrounds 


n/a 


TO 


Lapenta et al. (2013) 


21 


20 


A, C, S 


C3 and surrounds 


n/a 


TO 



General notes: DB= double-blind experimental design with regards to tDCS condition (not any included pharmacological arm); T(x)=time (in minutes) post stimulation that 
measurements were reported and pooled for this analysis. A=anode: C=cathode; S = sham. 

VEP notes: Bocci et al. (2013) reported conflicting information regarding stimulation strength and method of variability reporting. For pooling, we assumed and 1.5 niA 
current strength and variability reported as SEM. 

SEP notes: Hansen et al. (2011) was excluded from analysis due to reporting SDs over 100%. 



4. Discussion 

4.1. TMS discussion 

Of the six included TMS outcome measures, tDCS appears to 
modulate only MEP amplitude (although a significant difference 
was found between post-stimulation SICI scores for the cathode 
and anode conditions, neither of these scores differed significantly 
from sham nor from their own pre-stimulation levels). This finding 
is somewhat surprising because MEP amplitude is thought to re- 
flect corticocortical circuit excitability, which should also influence 
cSP, SICI, and ICE. Indeed, in the drug literature, any given inter- 
vention will rarely modulate MEP amplitude in isolation without 
effecting other TMS measure (for review: Paulus et al., 2008). A 
comparison of temporal cohorts suggests MEP modulation has 
been significantly decreasing since 2000. We will further discuss 
MEP amplitude modulation, the utility of this measure, and what 



the suggested temporal decrement might mean in the general 
discussion. Taken together, the evidence does not support the 
assertion that tDCS has a consistent effect on varied TMS outcome 
measures beyond MEP amplitude. 

4.2. ERP discussion 

tDCS does not appear to modulate any of the 6 pooled ERP 
outcome measures. Additionally, although no discernable numer- 
ical data was available for pooling, the written descriptions within 
included papers suggest tDCS does not modulate an additional 
9 ERP measures (PlOO VEP latency: low contrast stimuli, PlOO VEP 
latency: high contrast stimuli, N30 SEP amplitude, N30 SEP la- 
tency, Nl LEP Amplitude, Nl LEP latency, N2 LEP latency, P2 LEP 
latency, and LEP threshold). Taken together, the evidence does not 
support the assertion that tDCS has a consistent effect on varied 
ERP components. 






90 




70 


iselini 


50 


ta 


30 




10 




0 


3 


-10 


"5. 
E 


-30 


< 




SEP 


-50 


o 




z 


-70 




-90 



—Cathode 

—Anode 

-•-Sham 



Minutes Post Stimulation 



Minutes Post Stimulation 



Minutes Post Stimulation 



54 


D 


55 


in 


56 


n 


57 


I m 
1 «> 


58 


•S 20 


59 


* 0 

OP 


60 


1 -20 


61 




62 


z 

-80 


63 


-100 



■ Cathode 

■ Anode 

■ Sham 





Fig. 4. Effects of tDCS on various ERP amplitude measures. (A and B) Low-density, long duration ( > 7 min) tDCS exploring PlOO visual evoked potential amplitude to low- 
and high-contrast stimuli, respectively; (C) high-density, long duration tDCS exploring N20 somatosensory evoked potential amplitude; (d and E) low-density, long duration 
tDCS exploring N2 and P2 laser evoked potential amplitude, respectively; (F) long duration tDCS exploring Mu/Alpha ERD. Error bars: + 1 SD. 
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Table 14 

Numerical values for low-density, long duration tDCS on PlOO VEP amplitude (low 
and high contrast stimuli). 

iV Mean (% SD 95% u confidence 
of base- interval 
line) 

PlOO VEP amplitude: low and high contrast stimuli 
Low contrast 



Table 17 

Numerical values for long duration tDCS on Mu/Alpha ERD during motor imagery. 

IV Mean (% of baseline) SD 95% u confidence interval 



High contrast TO 



TO 


Cathode 


26 


14.77 


66.40 


-12.056 <u<41.596 




Anode 


26 


-8.62 


46.40 


-27.366<u< 10.126 




Sham 


- 


- 


- 


- 


TIO 


Cathode 


26 


17.85 


60.81 


-6.717 <u<42.417 




Anode 


26 


-5.62 


29.93 


- 17717 < u < 6.467 




Sham 


- 


- 


- 


- 


T20 


Cathode 


26 


-0.92 


50.15 


-21.181 <u< 19.341 




Anode 


26 


-9.77 


29.04 


-21.502 <u< 1.962 




Sham 










T30 


Cathode 


26 


-6.15 


31.43 


- 18.848 < u < 6.548 




Anode 


26 


-2.46 


21.85 


- 11.287 < u < 6.367 




Sham 










TO 


Cathode 


26 


13.30 


5782 


- 10.059 < u < 36.659 




Anode 


37 


1.63 


41.81 


-11.845 <u< 15.105 




Sham 










TIO 


Cathode 


26 


0 


41.91 


-16.932 <u< 16.932 




Anode 


26 


-9.08 


38.28 


-24.545 < IK 6.385 




Sham 










T20 


Cathode 


26 


-2.92 


35.19 


- 17137 < u < 11.297 




Anode 


26 


-3.69 


34.26 


- 17531 <u< 10.151 




Sham 










T30 


Cathode 


26 


-13 


34.73 


-27.031 <u< 1.031 




Anode 


26 


-5.39 


31.60 


- 18.156 < u < 7376 




Sham 











'Difference between either stimulation condition and sham. 
''Difference between stimulation conditions. 



Table 15 

Numerical values for high-density, long duration tDCS on N20 SEP amplitude. 



N Mean (% of baseline) 


SD 


95% u confidence interval 


N20 SEP amplitude 






TO Cathode 30 -5.33 


22.96 


-13.546 <u<2.886 


Anode 30 4.67 


25.26 


-4.369<u< 13.709 


Sham 10 8 


18.97 


-5.557 <u<21.557 


T20 Cathode 30 -6.17 


17.95 


-12.593<u<0.253 


Anode 30 1.33 


19.29 


-5.573<u<8.233 


Sham 10 0 


22.14 


-15.823 <u< 15.823 



'Difference between either stimulation condition and sham. 
''Difference between stimulation conditions. 

4.3. EEC power spectrum discussion 

tDCS does not appear to modulate any of the 4 pooled EEC 
power spectrum measures. It is possible different methodological 
parameters utilized in each study may have influenced this find- 
ing. However, as each included paper reported a significant finding 



Table 16 

Numerical values for low-density, long duration tDCS on N2 and P2 LEP amplitude. 



MU/alpha ERD: motor imagery 

TO Cathode 27 6.48 
Anode 35 6.43 
Sham 35 -10.26 



96.23 -31.670 <u<44.630 
83.65 -21.283 <u<34.143 
74.37 - 34.899 < u < 14.379 



■"Difference between either stimulation condition and sham. 
''Difference between stimulation conditions. 



on their respective outcome measure, one would expect to see an 
increased chance of finding a consistent/reliable effect. Taken to- 
gether, the evidence does not support the assertion that tDCS has a 
consistent effect on varied EEC oscillatoiy power measures; 
though it is important to emphasize that this conclusion is drawn 
solely from a qualitative comparison of similar studies. 

4.4. fMRl discussion „ 

toes does not appear to modulate voxel activation at Ml un- 
derneath the active electrode. In addition, although tDCS does not 
appear to modulate voxel activation at SMA following short- 
duration stimulation (either polarity) or cathodal long-duration 
stimulation, it is difficult to conclude the effect of anodal long- 
duration stimulation on this measure (as 2 studies reported a si- 
milar finding and only 1 reported the opposing finding). Taken 
together, however, the evidence does not support the assertion 
that tDCS has a consistent effect on cerebral blood flow measured 
by fMRl; though it is important to emphasize that this conclusion 
is drawn solely from a qualitative comparison of similar studies. 

4.5. General discussion 

We quantitatively and/or qualitatively reviewed the effects of 
transcranial direct current stimulation on a total of 30 neurophy- 
siologic measures across 4 different modalities. Our analysis of the 
literature shows that tDCS generates a significant and reliable ef- 
fect on only one measure: MEP amplitude (Table 22). The re- 
maining 29 measures did not reach significance or, in the case of 
those not amenable to pooling, displayed null or inconsistent 
results. 

One reason why many of these measures did not reach sig- 
nificance may be the lack of comparable research extant in the 
literature. It is certainly possible that with increased data (akin to 
the amount reporting MEP amplitude modulation), a number of 
additional physiologic outcome measures will prove amenable to 
tDCS modulation. Future replication research is certainly en- 
couraged to determine if this is the case. 







JV 


Mean (% of baseline) 


SD 


95% u confidence interval 


N2 and P2 LEP amplitude 












N2 amplitude TO 


Cathode 


32 


-29.06 


47.33 


-45.459 < IK -12.661 




Anode 


20 


-15.5 


43.25 


-35.712 <u<4.712 




Sham 


27 


-18.89 


47.67 


-37789 <u<0.009 


P2 amplitude TO 


Cathode 


32 


-19.94 


47.08 


-36.252 <u< -3.628 




Anode 


20 


-14.5 


37.80 


-32.203 <u<3.203 




Sham 


27 


-4.96 


34.27 


-18.546<u<8.626 



'Difference between either stimulation condition and sham. 
''Difference between stimulation conditions. 
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Table 18 

Studies exploring tlie effects tDCS on varied electroencepfialograpfiic (EEC) power spectrum measures. 



67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
91 
92 
93 
94 
95 
96 
97 
98 
99 
100 
101 
102 
103 
104 
105 
106 
107 
108 
109 
110 
111 
112 
113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 
131 
132 



Paper/* of studies 



N 



tDCS duration (min) 



tDCS polarities 



DB 



Active/reference 



Density (mA/cm^) 



Oscillation frequency power - at rest 

Ardolino et al. (2005) 6 C/5 S 10 

Matsumoto et al. (2010) 6 10 

Pellicciari et al. (2013) 16 13 

Notturno et al. (2013) 10 20 



C, S 
A, C, S 
A, C 
A, C, S 



Oscillation frequency power (and P300 ERP) - during N-back working memory task 



Keeser et al. (2011) 10 
Zaehle et al. (2011) 16 
Oscillation frequency characteristics - 
Wirth et al. (2011) 20 
Maeoka et al. (2012/x2) 15 



20 
15 

during picture tasks 
37 
20 



A, S 
A,C 

A,S 
A,S 



RMl/L orbit 


0.0429 


LMl/R orbit 


0.0286 


L Ml/R orbit 


0.04 


LMl/R orbit 


0.0286 


1 dlpfc/r orbit 


0.0571 


1 dlpfc/1 mastoid 


0.0286 


1 dlpfc/r shoulder 


0.0429 


1 dlpfc/r orbit 


0.0286 



General notes: DB= double-blind experimental design with regards to tDCS condition (not any included pharmacological arm); A^anode; C=cathode; S=sham; 

Ml ^primary motor cortex; DLPFC^dorolateral prefrontal cortex. 



Table 19 

Summary of findings from studies exploring the effects tDCS on varied electro- 
encephalographic (EEC) power spectrum measures. 

Gamma Beta Alpha Theta Delta 



Q40scillation frequency power - at rest 



t 



Cathode Ardolino et al. (2005) 
Matsumoto et al. 
(2010) 

Pelliciari et al. (2013) 
Notturno et al. (2013) 
Anode Matsumoto et al. 
(2010) 

Pelliciari et al. (2013) 

Notturno etal. (2013) - _ _ _ _ 

Oscillation frequency power - during N-back working memory task 
Anode Keeser et al. (2011 ) - - - - t 

Zaehle et al. (2011) - _ _ _ _ 

Oscillation frequency characteristics - during picture tasks 
Anode Wirth et al. (2011) - _ _ _ _ 

Maeoka et al. (2012) - _ _ _ _ 

i t - - 

Note; dash symbol represents no significant change. 

The complete power spectrum was measured for each included study. 



In addition, although MEP amplitude was the only outcome 
measure to achieve significance, there are a number of factors 
suggesting this result should be considered with caution. To begin, 
the majority of MEP amplitude studies in this review appear to 
have utilized TMS parameters that are not recognized to elicit 
reliable MEPs. Studies have long demonstrated large inter- and 
intra-subject variability of MEP amplitude; even for consecutively 
presented stimuli (Amassian et al., 1989; Kiers et al., 1993; Nielsen, 
1996; Roth and Magistris, 2008). To enhance the reliability of this 
measure, practitioners have argued that MEPs should be generated 
using a minimum TMS intensity of 130% rMT with a 6-10 s pause 
between each pulse (Brasil-Neto et al., 1992; Kiers et al., 1993). A 
look at Table 1 reveals that at least 71 included studies (out of 79) 
utilized a TMS intensity below the recommended 130% rMT (0.5- 
1 mV MEP amplitude equates to -115-125% rMT) and at least 50 
included studies delivered TMS pulses quicker than the re- 
commended 6-10 s ISI (0.25 Hz=one pulse every 4 sec). Further- 
more, it is worth noting human and animal research suggests that 
using 'stronger' TMS pulses increases the number and complexity 
of generated corticospinal volleys (Di Lazzaro et al., 2004; Ziemann 
and Rothwell, 2000; Amassian and Cracco, 1987). Pharmacological 
research suggests that later corticospinal volleys (not generated by 
lower intensity TMS) are more easily effected by and reflective of 
cortical excitability (specifically glutamatergic and GABAergic 
processes; for review: Paulus et al., 2008; also see Table S2 for a 



list of tDCS studies exploring the effects of tDCS on neuro- 
transmitter concentration). In fact, in a recent study, Lang et al. 
(2011) reported that tDCS effects are largest and most apparent in 
these later corticospinal volleys. Accordingly, the use of weaker 
TMS pulses to probe MEP amplitude modulation following tDCS 
may not only hinder reliability, but also may not completely probe 
the proposed mechanisms of change following stimulation. 

Additionally, as noted above, the effects of tDCS on MEP am- 
plitude have been significantly diminishing over recent years (most 
notably for anodal stimulation). Technological (both tDCS and 
TMS) and methodological advances since 2000 are likely to have 
reduced levels of noise which should, in turn, increase the relia- 
bility of outcome measures. The fact that MEP amplitude effects 
have been diminishing with these advances suggests that the ef- 
ficacy of tDCS may not be the sole explanation for changes in MEPs 
in the studies included in this review. 

It is possible that varying stimulation durations could account for 
the decrease in MEP amplitude effects since 2000. On average, more 
recent long-duration protocols have utilized a stimulation duration 
approximately 5 min longer than earlier studies {2000- 
2004= 10.50 min; 2005-2009= 10.62 min; 20W-Present= 15.32 min). 
According to the meta-plastic homeostatic theoiy (see Siebner et al., 
2004; Fricke et al., 2011), longer-duration stimulation should, in fact, 
reduce the modulatory response. However, our analysis also revealed a 
decrease in response for both veiy-shoit and short-duration protocols. 
Therefore, duration effects are unable to account for the complete 
diminution pattern we found. FurtheiTnore, as the number of studies 
which included a double-blind protocol and/or utilized neuronaviga- 
tion to guide TMS does not markedly vary between temporal cohorts, 
these paradigmatic consideration are also unable to account for this 
pattern. 

In a recent review paper, this research group pointed out five- 
major issues related to current tDCS research: namely, inter-sub- 
ject variability, intra-subject reliability, sham/blinding conditions, 
motor/cognitive interference, and electric current influences 
(Horvath et al., 2014). Each of these issues may provide potential 
explanations for the results obtained here. Regarding inter-subject 
variability, it is possible that many tDCS effects are unique to the 
individual being stimulated and, accordingly, these effects have 
been diminished or lost in pooled data. If this is the case, future 
experimentation exploring individualized response may better 
elucidate the efficacy and reliability of tDCS than a systematic 
review, such as this one. In addition, lack of accounting for hair 
thickness (a non-conductive barrier between tDCS and the head), 
individual neural architecture, and current flow dynamics may 
greatly influence tDCS effects which, in turn, may have influenced 
our findings. As noted above, some groups are now beginning to 
account for these effects. It is likely that the future adaptation of 
tDCS protocols to each individual may lead to a more 
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Table 20 

Studies exploring the effects tDCS on varied functional magnetic resonance imaging (fMRI) measures. 



Paper (# of studies) N tDCS duration tDCS Polarities DB Active/ 

reference 



Current density (niA/cm ) Task 



Compare 



Analysis notes 
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ROl: Ml activation/spread - 
Baudewig et al. 6 

(2001) 
Antal et al. (2011) 13 

ROl: Ml activation/spread - 
Stagg et al. (2009b) 15 

Jang et al. (2009) 7 

Saiote et al. (2013) 30 

ROl: SMA activation/spread 
Baudewig et al. 6 

(2001) 
Antal et al. (2011) 13 

Kwon and Jang (2011) 12 

ROl: SMA activation/spread 
Stagg et al. (2009b) 15 

Jang et al. (2009) 7 

Saiote et al. (2013) 30 



- during active motor task: short-duration 

5 min A, C - LMl/R orbit 



20 s 



A, C, control - LMl/R orbit 



- during active motor task: long-duration 

10 min A, C, S - LMl/R orbit 



20 min A, S 

10 min A, C, S 



LMl/R orbit 
LMl/R orbit 



— during active motor task: short-duration 

5 min A, C - LMl/R orbit 

20 s A, C. control - LMl/R orbit 

2 min A, S - IMl/r orbit 

- during active motor task: long-duration 

10 min A, C, S - LMl/R orbit 



20 min A, S 

10 min A, C. S 



LMl/R Orbit 
LMl/R orbit 



0.04 
0.0286 

0.0286 
0.0286 
0.0286 

0.04 

0.0286 

0.0286 

0.0286 
0.0286 
0.0286 



Finger tapping 
Finger tapping 



Sequential finger 
press 

Grasp/release 
Fist squeezing 



Finger tapping 
Finger tapping 
Grasp/release 



Sequential finger 
press 

Grasp/release 
Fist squeezing 



Pre- to post-stim 

Post-stim active vs 
control 

Post-stim active vs 
sham 

Post-stim active vs 
sham 

Post-stim active vs 
sham 

Pre- to post-stim 

Post-stim active vs 
control 

Post-stim active vs 
sham 

Post-stim active vs 
sham 

Post-stim active vs 
sham 

Post-stim active vs 
sham 



r:COA=p < 0.001 Spread=p < 0.05 (software: In-house) 

vol RFX: corrected threshold=(;(FDR)=0.05(software: 
BrainVoyagerQX) 

ROl MFX: corrected threshold =Z> 2.0, p < 0.01(software: 
FMRIB) 

ROl RFX: uncorrected threshold =p < 0.001 (software: SPM2) 
MFX: corrected threshold =Z< 2.3. p < 0.05(software: FEAT) 



r:COA=p < 0.001 spread=p < 0.05 (software: In-house) 

vol RFX: corrected threshold = (;(FDR)=0.05(software: 
BrainVoyagerQX) 

vol RFX: uncorrected Threshold =p < 0.001 (software: SPM8) 



ROl MFX: corrected Threshold =Z > 2.0, p < 0.01 (software: 
FMRIB) 

ROl RFX:Uncorrected Threshold =p < 0.001 (software: SPM2) 
MFX: corrected threshold =Z< 2.3, p < 0.05(software: FEAT) 



3 



General notes: DB= double-blind experimental design with regards to tDCS condition (not any included pharmacological arm); R01=region of interest: A^anode; C^cathode; S=sham; Ml = primary motor cortex; r=correlation 
coefficient; COA=center of activation; IND^individual noise distribution; V01=volume of interest; RFX=random effects analysis; qFDR^false discovery rate q-value; MFX = mixed-effects modelROl: Ml activation/spread: during 5 liH 
active motor task. 
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Table 21 

Summary of findings from studies exploring the effects tDCS on varied functional 
magnetic resonance imaging (fMRl) measures. 



Peak voxel 
activation 



Cluster size (# of 
voxels) 



Table 22 

Summary of the effect of tDCS on each neurophysiological outcome measure ex- 
plored in this systematic review. 

Numerically analyzable outcome measures 



QI3MI activation during active motor task: short-duration 

Cathode Baudewig et al. (2001 ) - 

Antal et al. (2011) 
Anode Baudewig et al. (2001) - 

Antal et al. (2011) 
Ml activation — during active motor task: long-duration 
Cathode Stagg et al. (2009a, 

2009b) 

Saiote et al. (2013) 
Anode Jang et al. (2009) 

Stagg et al. (2009) t T 

Saiote et al. (2013) 
SMA activation — during active motor task: short-duration 
Cathode Baudewig et al. (2001 ) i I 

Antal et al. (2011) 
Anode Baudewig et al. (2001) 

Antal et al. (2011) I 4, 

Kwon and Jang (2011) t t 
SMA activation — during active motor task: long-duration 
Cathode Stagg et al. (2009a, 

2009b) 

Saiote et al. (2013) 
Anode Jang et al. (2009) t T 

Stagg et al. (2009a, t T 
2009b) 

Saiote et al. (2013) 



Note: dash symbol represents no significant change. 



comprehensive and consistent view of tDCS effects. 

Finally, the lacl< of sham-controlled, double-blind, and neuio- 
navigation-guided (NN) TMS studies is somewliat alarming. Of 106 
TMS Studies included in the analysis, only 18 included a sham/ 
control condition, only 7 utilized a double-blind study design, and 
only 5 reported using NN to guide TMS delivery. Although the lack 
of studies using NN may be understandable (as the technology is 
rather novel), it is worth noting that growing body of evidence 
demonstrates that NN allows for millimeter precision whilst 
"standard" TMS positioning can engender significant inter-pulse 
coil movement (see Gugino et al., 2001 ; Schonfeldt-Lecuona et al., 
2005; Sparing et al., 2008; Sack et al., 2009; ) In fact, Ahdab et al. 
(2010) have argued that the inaccuracy of non-NN TMS is so great 
that results obtained using "standard" positioning techniques 
should be interpreted with caution. The lack of studies utilizing a 
double-blind paradigm may also be somewhat understandable, as 
research has demonstrated it is quite difficult to blind the tDCS 
practitioner, especially if he/she is the same person delivering TMS 
or setting up EEC electrodes (a necessity in many smaller research 
groups: O'Connell et al., 2012). However, the lack of studies in- 
cluding a sham condition is harder to understand as this is clearly 
established scientific best-practice. Future research would benefit 
from the inclusion of a sham condition, blinding of both the par- 
ticipant and experimenter to the tDCS condition being utilized, 
and employment of modern frameless stereotactic NN systems 
(when utilizing TMS) to ensure all results are accurate and 
unbiased. 

4.6. Limitations 

A major limitation of this systematic review is the lack of 
comparable research extant in the literature. Beyond the many 
neurophysiological outcome measures in the literature not ex- 
plored by at least two different research groups (see Table SI in 
Supplementary material), many replicated outcome measures 



tDCS outcome measure (modality) 

MEP amplitude (TMS) 

Resting motor threshold (TMS) 

Active motor threshold (TMS) 

Cortical silent period (TMS) 

Short interval cortical inhibition (TMS) 

Intracortical facilitation (TMS) 

PlOO visual evoked potential amplitude - low contrast 
stimuli (ERP) 

PlOO visual evoked potential amplitude - high contrast 
stimuli (ERP) 

N20 somatosensory evoked potential amplitude (ERP) 

N2 laser evoked potential amplitude (ERP) 

P2 laser evoked potential amplitude (ERP) 

MU event related displacement (ERP/EEG power spectrum) 

Non-numerically analyzable outcome measures 

tDCS outcome measure 

PlOO visual evoked potential latency - low contrast stimuli 
(ERP) 

PlOO visual evoked potential latency - high contrast stimuli 
(ERP) 

N30 somatosensory evoked potential amplitude (ERP) 

N30 somatosensory evoked potential latency (ERP) 

Nl laser evoked potential amplitude (ERP) 

Nl laser evoked potential latency (ERP) 

N2 laser evoked potential latency (ERP) 

P2 Laser Evoked Potential Latency (ERP) 

Laser evoked potential threshold (ERP) 

Oscillatory frequency power at rest (EEC power spectrum) 

Oscillatory frequency power during N-back task (EEC 

power spectrum) 
P300 potential amplitude during N-back task (ERP/EEG 

power spectrum) 
P300 potential latency during N-back task (ERP/EEG power 

spectrum) 

Oscillatory frequency power during picture tasks (EEG 

power spectrum) 
Ml activation during active motor task - short duration 

stimulation (fMRl) 
Ml activation during active motor task - long duration 

stimulation (fMRl) 
SMA activation during active motor task - short duration 

stimulation (fMRl) 
SMA activation during active motor task - long duration 

stimulation (fMRl) 



Significant effect 
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have utilized different values of the four adjustable tDCS para- 
meters (electrode location, current density, stimulation duration, 
and stimulation-to-task relationship). Due to this variability, there 
simply is a dearth of meaningfully comparable tDCS studies. It is 
certainly possible that, with increased data, a number of additional 
measures (beyond MEP amplitude) will reach significance. As this 
field evolves, it may be beneficial to consider standardized pro- 
tocols and attempt to directly replicate work already in the pub- 
lished literature. 

Another limitation concerns the inability to derive a mean- 
ingful effect size from the neurophysiological data. As noted in the 
introduction, the majority of neurophysiological studies included 
in this review did not use a control or sham condition; rather, 
many simply compared either varied active stimulation protocols 
or pre- and post-active stimulation measures (for comment: Hor- 
vath et al., 2014). Due to an inability to derive meaningful effect 
sizes, we were unable to conduct a true meta-analysis. Accord- 
ingly, the quantitative analyses used in this systematic review 
(derived from weighted means and pooled standard deviations) 
must be interpreted with caution. It is certainly possible that, had 
it been feasible to derive meaningful effect sizes, a number of 
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measures may have reached significance. Again, as this field 
evolves, it may be beneficial to include effective control conditions 
in eveiy study so that the effects of stimulation can be differ- 
entiated from natural neurophysiological fluctuation and a more 
accurate effect of stimulation can be determined. 

Along similar lines, a number of papers included in this sys- 
tematic review were not amenable to pooling due to no hard-data 
being apparent in the final publication/s (either in the form of 
numerical tables or graphs). It is possible, had all researchers in- 
cluded hard-numbers (even for measures that did not reach sig- 
nificance in the original study), a number of outcome measures 
may have reached significance in this review. Accordingly, future 
researchers should be encouraged to include discernable data so 
that future analyses can include more data and develop a more 
accurate picture of tDCS effect. 

Another limitation of this study is its exclusion of data obtained 
in non-healthy adult populations. There is a large body of work 
exploring the neurophysiologic effects in clinical populations (for 
review: Brunoni et al., 2012), elderly population (e.g. - Berryhill 
and Jones, 2012; Kar and Wright, 2014), and animal models (for 
review: Brunoni et al., 2012). Accordingly, the results and con- 
clusions reported in this paper cannot and should not be extra- 
polated to these other groups. 

Finally, there is a rather robust literature suggesting tDCS can 
modulate cognitions and performance on a number of behavioral 
tasks. Although the lack of significant neurological effects found in this 
review may raise questions concerning the mechanisms by which 
tDCS modulates behavior, it does not speak directly (positively nor 
negatively) to behavioral modulation itself As such, it would be a 
mistake to extend the findings of this analysis to the whole of the tDCS 
literature; rather, similar, robust systematic reviews of the cognitive 
and behavioral data are necessaiy before far reaching comments about 
this device can be made. This group has undertaken such reviews and 
the manuscripts are currently in preparation. 



5. Conclusion 

In the end, tDCS does not appear to generate a reliable neuro- 
physiologic effect. Although MEP amplitude does appear to be sensi- 
tive to tDCS modulation, this effect has been significantly decreasing 
since 2000 and other, more reliable TMS measures believed to rely on 
similar neural mechanisms (e.g. SICI, ICF, and cSP) have all shown no 
effect to toes. Tliese findings raise questions about the mechanistic 
framework of tDCS. If future research continues to generate evidence 
consistent with the conclusion that tDCS generates little-to-no reliable 
neurophysiological effect, it will be important to begin reconsidering 
the current mechanistic foundations used to explain the effects of 
tDCS. It is hoped that greater consistency in the utilization of stimu- 
lation protocols, the use of sham conditions, and the statistical re- 
porting of all outcome measures assessed (regardless of significance) 
will allow for a more comprehensive assessment of the utility of tDCS 
in the future. 



Uncited references 

Cogiamanian et al. (2007), Kidgell et al. (2013b), Madhavan and 
Stinear (2010) and Rosier et al. (2008). 



Appendix A. Supplementary material 

Supplementary data associated with this article can be found in 
the online version at http://dx.doi.0rg/lO.lOl6/j.neuropsychologia. 
2014.11.021. 
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